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One  of  the  essential  problems  in  the  design  of 
mechanisms  is  the  force  and  torque  transmission.  The  trans- 
mission problem  for  both  planar  mechanisms  and  the  spatial 
mechanisms  containing  lower  pairs  has  been  discussed  exten- 
sively. However,  this  is  not  so  for  the  spatial  mechanisms 
containing  higher  pairs.  The  objective  of  this  work  is  to 
study  the  transmission  problem  of  the  particular  mechanisms. 

The  transmission  characteristics  of  spatial  mechanisms 
containing  higher  pair  joints  are  investigated  by  way  of  the 
deviation  angles  and  the  transmissivities,  which  are  ex- 
pressed by  way  of  geometric  parameters.  Sliding  friction 
forces  are  included  in  the  force  and  torque  analysis  due  to 
their  significant  affect  on  the  direct  contact  between  the 
higher  pair  elements.  Thus,  mechanical  advantage  and  in- 
stantaneous efficiency  are  also  investigated  for  obtaining 
the  efficiency  of  the  mechanisms.  In  order  to  maintain 
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continuous  motion,  sliding  friction  forces  must  be  overcome 
by  the  effective  force  produced  by  the  input  power.  Thus,  a 
geometric  relationship  between  the  friction  force,  which  is 
opposes  the  direction  of  the  relative  velocity  on  the  con- 
tact plane,  and  the  effective  force  on  the  same  contact 
plane  must  be  satisfied  in  addition  to  the  deviation  angles. 

Higher  pair  elements  may  have  point  or  line  contact. 

To  reduce  the  complexity  of  the  formulations  caused  by  dis- 
tributed friction  forces,  line  contact  is  modeled  by  two 
concentrated  normal  force  applied  at  the  ends  of  the  contact 
line.  Similarly,  for  the  plane  pair,  a lower  pair,  which  is 
considered  as  a higher  pair  despite  of  the  surface  contact, 
is  modeled  by  three  concentrated  normal  forces  acting  on 
three  of  the  corners  of  the  contact  plane. 

Numerical  examples  are  given  for  the  design  of 
mechanisms  containing  one-point,  line,  and  surface  contact 
higher  pairs,  respectively,  which  includes  the  consideration 
of  optimum  deviation  angles  and  minimized  peak  acceleration 
to  improved  the  high  speed  performance. 
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CHAPTER  1 
INTRODUCTION 

1 . 1 Motivation 

The  kinematic  synthesis  and  analysis  of  closed  loop 
spatial  mechanisms  have  been  studied  extensively  on 
mechanisms  with  lower  pairs,  but  less  so  on  mechanisms  con- 
taining higher  pairs.  The  use  of  higher  pairs  leads  to  more 
compact  spatial  mechanisms,  which  in  general  can  be  designed 
to  produce  better  motion  approximations  and  to  have  better 
dynamic  performance  than  lower-pair  mechanisms. 

One  of  the  essential  problems  in  the  design  of 
mechanisms  is  the  force  and  torque  transmission  which  is 
also  a criterion  to  avoid  the  "lock  up"  positions 
(Sutherland,  1981)  or  uncertainty  configurations  as  defined 
by  Hunt  (1978).  The  transmission  problem  for  both  planar 
mechanisms  and  spatial  mechanisms  containing  lower  pairs  has 
been  discussed  extensively  in  existing  literature.  However, 
this  is  not  the  case  for  spatial  mechanisms  containing 
higher  pairs.  Therefore,  the  objective  of  this  work  is  to 
study  the  transmission  problem  of  spatial  mechanisms  with 
higher  pairs. 

Instead  of  using  transmission  angles  as  defined  in  the 
design  of  planar  mechanisms  and  some  spherical  and  also  in 
spatial  mechanisms  containing  lower  pairs  only,  deviation 
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angles  are  used  as  the  design  criteria.  The  deviation 
angles  are  defined  as  the  angles  between  the  input  force  and 
the  resistant  force  at  each  joint  in  the  sequence  of  force 
transmission  which  describe  the  transmitted  force  at  the 
floating  joints.  The  deviation  angles  for  moment  transmis- 
sion are  also  defined  similarly  to  those  for  the  force 
transmission.  For  each  floating  joint  there  will  be  devia- 
tion angles  for  the  transmission  of  internal  force  and  mo- 
ment, i.e.,  the  force  and  moment  transmitted  between  the 
pair  elements.  If  the  joint  connects  the  coupler  link  and 
the  output  link,  then  there  will  also  be  deviation  angles 
for  the  force  and  moment  to  produce  the  output  force  and 
output  moment.  The  force  and  moment  transmissivities  are 
introduced  to  take  all  deviation  angles  in  sequence  into  ac- 
count. For  three-link  mechanisms,  the  transmissivities  are 
classified  as  force-to-f orce  and  f orce-to-torque . These  are 
the  indexes  of  how  well  the  coupler  force  produces  the  out- 
put translational  and  rotational  motions,  respectively.  The 
torque-to-torque  transmissivity  is  the  index  of  how  well  the 
coupler  produces  the  rotational  motion  of  the  output  joint. 
For  four-link  mechanisms,  the  coupler  force  will  produce  a 
moment  at  the  joint  connecting  the  coupler  link  and  the  out- 
put link.  Thus  a f orce-torque-to-torque  transmissivity  is 
needed  to  show  how  this  moment  produces  the  output  rota- 
tional motion.  To  take  the  frictional  effects  into  account, 
the  mechanical  advantage  and  the  instantaneous  efficiency 
are  also  investigated.  The  mechanical  advantage  is  defined 
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as  the  ratio  between  the  output  force  exerted  by  the  output 
link  and  the  input  force  applied  to  input  link.  The  instan- 
taneous efficiency  is  defined  as  the  ratio  between  the  out- 
put power  and  the  input  power.  When  there  is  no  energy  loss 
due  to  the  frictional  effects,  the  instantaneous  efficiency 
should  be  unity.  The  mechanical  advantage  and  the  instan- 
taneous efficiency  of  the  mechanism  are  affected  by  the 
power  loss  due  to  frictional  effects  between  pair  elements. 
In  the  consideration  of  high  speed  performance,  the  dynamic 
behavior  of  the  mechanism  should  also  be  of  concern.  This 
can  be  evaluated  by  the  angular  and  linear  accelerations  of 
the  mechanism  links.  The  output  acceleration  ratios  can  be 
defined  as  the  ratio  between  the  maximum  and  minimum  output 
angular  and  linear  accelerations,  which  depend  on  the  type 
of  the  output  joint.  If  the  output  link  has  only  angular 
motion,  the  acceleration  ratio  may  be  defined  as  the  ratio 
between  the  output  angular  acceleration  and  the  square  of 
the  constant  input  angular  velocity,  as  used  in  tandem 
mechanisms  (Hain,  1971).  Since  the  input  angular  velocity 
is  assumed  to  be  constant  in  this  work,  therefore,  the  ac- 
celeration ratio  will  be  replaced  by  the  magnitude  of  the 
output  angular  acceleration  itself.  It  also  may  be  desired 
to  minimize  the  deviation  between  the  maximum  and  minimum 
output  angular  or  linear  accelerations. 

Once  the  transmission  criteria  are  defined,  the  trans- 
mission characteristics  of  a certain  mechanism  can  be 
evaluated  by  way  of  the  deviation  angles , the  force  and  mo- 
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ment  transmissivities,  the  mechanical  advantage,  and  the  in- 
stantaneous efficiency.  The  following  step  then  is  to  op- 
timize the  mechanism  to  create  better  transmission  and  bet- 
ter dynamic  performance.  The  other  critical  factors  to  be 
considered  are 

1.  the  rotatablity  of  the  mechanism; 

2.  the  physical  size  of  the  higher  pair,  which  may 
create  a restriction  of  the  relative  motion  between 
two  pair  elements;  and 

3.  the  avoidance  of  possible  redundant  translational 
freedoms,  i.e.,  the  possible  collinearity  of  sliding 
axes . 


1 . 2 Literature  Review 

1.2.1  Kinematics  of  the  Mechanisms  Containing  Higher  Pairs 
Kohli  and  Soni  (1973a,  b)  use  successive  screw  dis- 
placements to  conduct  kinematic  analysis  of  single-loop  and 
two-loop  spatial  mechanisms  involving  lower  pairs.  Kohli 
and  Soni  (1975)  and  Singh  and  Kohli  (1979)  developed  pair 
geometry  constraint  equations  for  lower  pairs  such  as 
Revolute  (R),  Prismatic  (P),  Helical  (H),  Sphere  (S),  and 
Cylinder  (C)  pairs,  and  for  higher  pairs  such  as  Slotted 
Sphere  (SL),  Sphere-in-slotted  Cylinder  (Sc),  Sphere-Plane 
(Sp),  Sphere-Groove  (Sg)  and  Cam  pairs.  Figure  (1-1)  il- 
lustrates these  higher  pairs.  The  higher  pairs  can  also  be 
represented  by  various  combinations  of  lower  pairs  as  seen 
in  the  work  of  Chen  and  Roth  (1969).  Litvin  and  Gutman 
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Figure  1-1  Pair  Types. 
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Figure  1-1  continued. 
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(1980)  use  matrix  methods  for  a three-link  mechanism  with  an 
intermediate  higher  pair  of  two  crossed  cylinders.  Sandor, 
Kohli,  Hernandez  and  Ghosal  (1981)  and  then  Ghosal,  Kohli 
and  Sandor  (1982)  applied  higher  pair  geometry  constraint 
equations  to  analyze  three-link  spatial  mechanisms  that  con- 
tain Sp  and  Sg  pairs  and  four-link  mechanisms  that  contain 
Sg  and  Sc  pairs,  respectively.  A more  complete  kinematic 
analysis  is  done  by  Ghosal  in  his  thesis  (1982)  including 
three,  four  and  five  link  mechanisms.  Hernandez  presented 
the  kinematic  synthesis  and  analysis  of  three-link  spatial 
function  generators  containing  higher  pairs  including  pair 
element  design  (1983).  Miesch  (1984)  gave  optimal  design  of 
the  R-Sp-R  function  generator.  However,  transmission 
characteristics  have  not  been  discussed  fully  as  yet  for 
higher  pair  mechanisms,  although  this  is  imperative  in  the 
design  of  such  mechanisms. 

1.2.2  Coulomb  Friction  in  Higher-pair  Mechanisms 

The  contact  between  higher  pair  elements  is  likely  to 
be  influenced  by  friction  forces  that  are  often  distributed 
along  a line  or  curve,  rather  than  being  a concentrated 
force  acting  at  a single  point.  Although  the  friction  can 
be  neglected  in  the  case  of  well  lubricated  bearings  in 
lower  pairs,  for  some  circumstances  the  lubrication  is  poor 
or  not  allowed  at  all  (for  example,  in  space  applications); 
then  the  friction  is  to  be  taken  into  account.  There  are 
not  many  publications  in  mechanism  dynamics  that  consider 
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the  effects  of  friction  in  lower-pair  bearings.  Keler 
(1973)  derived  the  static  forces,  including  friction  of 
lower  and  some  higher  pairs  using  dual  numbers  based  on  an 
iterative  technique,  that  is  applicable  to  all  single-loop 
mechanisms  which  are  kinematically  and  statically  deter- 
minate. Sadler  (1973)  uses  the  concept  of  the  friction 
circle  to  investigate  the  effects  of  Coulomb  friction  in  a 
pin  joint  with  equilibrium  equations  in  the  dynamic  analysis 
of  four  bar  mechanisms.  Bagci  (1975a)  presented  the  dynamic 
motion  analysis  of  statically  determinate  planar  mechanisms 
via  the  joint-force  analysis  including  Coulomb  and  viscous 
damping.  The  methods  used  by  Sadler  (1973)  and  Bagci 
(1975a)  parallel  the  graphical  approach.  Imam,  Skreiner, 
and  Sadler  (1981)  utilize  similar  method  but  they  add  the 
approximate  closed-form  solution  to  force  analysis  in  the 
application  of  circuit-breaker  mechanisms.  Recently,  Bagci 
and  Khosravi  (1986)  extended  their  previous  work  (1975)  to 
the  application  to  spatial  mechanisms  with  lower  pairs,  such 
as  revolute,  cylinder,  prismatic,  and  sphere  pairs.  In  the 
above  work,  all  line  contacts  are  treated  as  contacts  at  a 
single  point  with  the  assumption  that  only  one  concentrated 
reaction  force  exists  at  the  center  of  gravity  of  the  con- 
tact line.  This  is  allowable  because  of  the  relative  linear 
translation  between  pair  elements,  and  does  not  involve  the 
relative  rotation  on  the  tangential  plane  of  the  contact. 
This  assumption  is  also  applicable  to  higher  pairs  with 
one-point  contact,  but  not  for  higher  pairs  with  more-than- 
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one-point  contact,  where  the  relative  velocity  is  different 
for  each  contact  point  along  the  line  of  contact. 

Consideration  of  Coulomb  friction  is  also  seen  in  the 
design  of  cam  systems.  Pisano  (1984)  modeled  a high-speed 
cam  system  using  Coulomb  friction  alone  at  the  rocker  arm 
pivot.  The  static  and  dynamic  coefficients  of  Coulomb  fric- 
tion are  calculated  directly  from  experimental  data  taken  at 
low  camshaft  speeds.  The  incorporation  of  Coulomb  friction 
into  the  model  gives  very  good  agreement  with  experiment. 
Chew  (1985)  used  the  concept  of  rolling-contact  to  reduce 
friction  at  highly  loaded  rocker-pivot  without  the  use  of 
rolling  element  bearings. 

Another  subject  is  the  friction  in  the  sliding  surface 
of  slider  mechanisms.  The  study  in  this  area  is  more 
dynamic,  which  may  include  surface  combination  (different 
sliding  surfaces  in  contact),  surface  finish,  viscosity  of 
the  lubricant,  the  composition  of  the  lubricant,  the 
geometry  of  the  lubricating  grooves  and  the  nominal  area  of 
contact  with  apparent  pressure  on  the  slideway  (Bell  & Bur- 
dekin,  1966).  Obviously,  these  subjects  are  outside  the 
scope  of  what  is  considered  in  this  work. 

Two  technigues  for  reducing  the  friction  between 
kinematic  pairs  which  connect  a movable  link  to  a frame  are 
described  by  Litvin  and  Coy  (1984).  One  way  is  arranging 
the  design  so  that  the  relative  velocities  of  the  link  with 
respect  to  two  grounded  kinematic  pairs  are  egual  in  mag- 
nitude and  opposite  in  direction.  The  other  way  is  moving 
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the  link  in  a direction  perpendicular  to  the  direction  of 
the  working  motion.  Although  these  techniques  are  not  ap- 
plicable in  higher  pairs,  since  higher  pairs  are  usually  at 
floating  joints,  it  may  be  applied  to  grounded  lower  pairs 
if  frictional  effects  are  considered. 

1.2.3  Coulomb  Friction  Coefficient 

The  static  coefficient  of  Coulomb  friction  in  the 
kinematic  pairs  of  a mechanism  is  usually  used  when  the 
mechanism  is  about  to  move.  When  the  mechanism  is  already 
in  motion,  the  kinetic  coefficient  of  friction  is  applied, 
which  in  general  is  about  75  percent  of  the  static  coeffi- 
cient (Bagci,  1975a). 

In  sliding  motion,  it  may  happen  that  the  magnitude  of 
the  relative  velocity  goes  to  zero,  for  instance  at  the  mo- 
ment when  the  mechanism  is  starting  up,  reversing,  or  a 
position  is  broken.  To  deal  with  zero  relative  velocities, 
several  models  have  been  seen  in  the  design  of  mechanisms 
containing  a sliding  joint.  The  most  straightforward  model 
is  that  the  direction  of  the  friction  force,  F,  depends  on 
the  sign  of  the  relative  velocity,  Vs,  which  is  expressed 

F = u|N| (-sign(Vs) ) (1-1) 

where  |N|  is  the  absolute  magnitude  of  the  normal  force  ap- 
plied at  the  contact  point  p.  Figure  (l-2a)  illustrates 
Equation  (1-1)  where  the  effect  of  zero  relative  velocity  is 
not  included.  Bell  and  Bardekins ' experimental  measurements 
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Figure  1 2 Approximations  of  static  friction  coefficient  at 
zero  relative  velocity. 
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(1966)  show  that  Equation  (1-1)  is  especially  good  for 
moderate  sliding  velocities  and  cyclic  conditions.  Some 
other  models  are  discussed  by  Rooney  and  Deravi  (1982), 
which  includes  a straight  line  approximation  when  the  rela- 
tive velocity  goes  to  zero,  as  shown  in  Figure  (l-2b),  and 
the  one  proposed  by  Rooney  and  Deravi,  which  has  a nonlinear 
approximation  nears  zero  velocity,  as  shown  in  Figure  (1- 
2c).  Keler  (1973)  suggests  using  the  relative  acceleration 
vector  as  the  direction  of  the  sliding  friction  force  and 
using  the  adhesion  coefficient  instead  of  the  frictional 
coefficient.  In  the  present  work.  Equation  (1-1)  is  ap- 
plied. However,  generally  the  sliding  velocity  on  the  con- 
tact plane  of  the  higher  pairs  does  not  go  to  zero  and 
changes  to  the  opposite  direction  during  the  motion  of  the 
mechanism,  i.e.,  generally  the  sliding  motion  is  not  along  a 
straight  line,  and  therefore  application  of  Equation  (1-1) 
is  not  quite  straightforward. 

1-2.4  Transmission  Criteria  for  Linkages 

For  planar  and  spherical  four-bar  mechanisms,  transmis- 
sion criteria  are  defined  by  the  well-known  transmission 
a^9ls , i.e. , the  angle  between  the  coupler  and  the  output 
link.  The  first  application  of  optimal  design  including  the 
transmission  angle  in  planar  mechanisms  probably  is  the  work 
done  by  Roth,  Freudenstein,  and  Sandor  (1962).  In  the  ensu- 
ing two  decades  and  more,  voluminous  publications  consider- 
ing the  optimal  transmission  angle  were  published  ( Shoup  & 
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Pelan,  1971;  Freudenstein  & Primrose,  1972;  Bagci,  1975b; 
Gupta,  1977  & 1978;  Freudenstein  & Chew,  1979;  and  Tsai, 
1983).  Most  of  these  works  use  the  transmission  angle  to 
define  the  transmission  characteristics,  which  is  generally 
adequate  for  planar  mechanisms.  Graphical  procedures  for 
the  optimization  of  the  transmission  angle  in  planar 
mechanisms  are  seen  in  Hall's  work  (1966)  and  also  in  the 
work  by  Sun  and  Waldron  (1981). 

The  transmission  angle  is  also  used  in  the  design  of 
spherical  mechanisms  (Dukkipati  & Soni,  1975;  Huang,  Pamidi, 
& Soni,  1970,  Soni  & Harrisberger , 1967;  Hossne,  1979; 
Chiang,  1986a,  b;  Huang  & Soni,  1969)  and  in  spatial 
mechanisms  (Soni  & Harrisberger,  1969,  1971).  Brasfield  and 
Bagci  (1975b)  claimed  that  the  use  of  transmission  angle 
defined  similarly  to  those  used  for  planar  mechanisms  may 
not  be  appropriate  for  spherical  and  spatial  mechanisms. 

To  express  the  transmission  criteria  for  spatial 
mechanisms  with  lower  pairs,  various  works  have  used  dif- 
ferent approaches.  Hamid  and  Soni  use  the  definition  of 
transmission  ratio,  which  is  the  ratio  of  the  force  exerted 
by  the  coupler  link  to  the  applied  force  received  by  the 
output  link.  They  applied  these  in  the  design  of  RSSR 
(Hamid  & Soni,  1971)  and  RSSP  (Hamid  & Soni,  1972) 
mechanisms.  In  the  publication  of  Shoup,  Steffen  and 
Weatherford  (1972),  two  deviation  angles  are  defined  as  the 
force  and  torque  effectiveness,  which  are  1)  the  angle  be- 
tween the  direction  of  the  force  exerted  by  the  coupler  on 
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the  output  link  and  the  direction  of  the  velocity  of  their 
joint  and  2)  the  angle  between  the  axis  of  the  torque  ex- 
erted by  the  coupler  on  the  output  link  and  the  direction  of 
the  axis  of  the  angular  velocity  of  the  output  link.  Bagci 
(1971)  uses  3x3  screw  matrices  to  investigate  the  force-to- 
force,  f orce-to-torque , torque-to-torque , and  force  torque- 
to-torque  transmissivities  which  are  solved  through  dual 
equilibrium  equations.  Transmission  criteria  are  thus 
determined  by  how  the  force  and  moment  exerted  by  the 
coupler  affect  the  output  force  and  moment.  Bagci  (1971) 
also  indicated  that  the  transmissivities  can  be  used  to 
avoid  passive  or  redundant  freedom.  This  is  because  the 
transmissivities  vanish  when  there  is  passive  or  redundant 
freedom.  The  torque-to-torque  transmission  ratio  concept  is 
also  used  in  the  later  work  of  Bagci  (1973)  for  the  spheri- 
cal crank-rocker  mechanism.  Also,  in  the  design  of  RCRRC 
mechanism,  the  transmissivities  are  applied  by  Lee  and  Bagci 
(1975)  again.  Sutherland  (1973)  established  a transmission 
index  by  using  the  definition  of  wrenches  and  screws  related 
to  the  mechanical  error  in  a linkage.  The  transmission  in- 
dex can  be  used  as  an  index  of  the  quality  of  motion 
(Sutherland,  1981).  Litvin  (1980)  uses  dual  vectors  to 
develop  the  force  transmission  criteria  for  spatial  linkages 
of  any  structure  where  he  defines  the  transmission  criteria 
by  the  ratios  of  the  force  and  moment  of  every  joint  to  the 
resisting  moment  at  the  joint.  A definition  of  transmission 
angle  similar  to  the  definition  of  transmission  ratio  is 
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used  by  Hamid  and  Soni  (1971,  1972),  and  it  is  also  seen  in 
the  works  of  Gupta  and  Kazerounian  (1983),  Rao  and 
Lakshminarayana  (1984a,  b,  c),  and  Soylemez  and  Freudenstein 
(1982).  In  these  works,  the  transmission  ratios  are  used  to 
denote  the  ratios  between  the  force  and  moment  exerted  by 
the  coupler  to  the  output  force  and  moment. 

The  use  of  transmission  angle  or  ratio  is  well  suited 
for  single  force  transmitting  mechanisms  such  as  the  RSSR 
mechanism.  For  mechanisms  containing  higher  pairs,  which 
may  have  both  force  and  torque  transmitted,  the  definition 
similar  to  those  given  by  Bagci  (1971)  or  Litvin  (1980) 
would  be  more  appropriate.  Similar  definition  of  transmis- 
sivities to  that  used  by  Bagci  (1971)  is  applied  with  a dif- 
ferent approach  in  this  work.  Funabashi  and  Freudenstein 
(1979a,  b)  gave  some  transmission  criteria  for  high  speed 
planar  and  spherical  crank-and-rocker  linkages.  These 
criteria  may  be  also  applicable  in  the  design  of  those  par- 
ticular mechanisms  discussed  in  this  work  which  are  usually 
used  for  high  speed  applications.  This  is  similar  to  the 
design  of  cam  mechanisms,  where  the  ratio  of  the  output  an- 
gular acceleration  to  the  square  of  the  input  angular 
velocity  is  considered  (Hain,  1971). 

1 • 3 Objectives  and  Summary 

Most  works  mentioned  previously  deal  with  lower  pairs 
without  Coulomb  friction  while  the  purpose  of  this  work  is 
to  determine  the  transmission  criteria  of  spatial  mechanisms 
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including  Coulomb-frictional  effects  in  higher  pair  joint. 
The  friction  forces  do  not  affect  the  deviation  angles,  but 
they  do  affect  the  reaction  forces  and  moments.  Thus,  be- 
side the  deviation  angles,  mechanical  advantage  and  instan- 
taneous efficiency  are  also  needed.  Developing  definitions 
of  better  transmission  criteria  will  be  of  benefit  in  the 
seeking  of  optimal  designs  of  spatial  mechanisms  with 
higher-pair  joints. 

In  this  chapter  we  have  reviewed  the  works  related  to 
the  topics  of  this  dissertation  and  have  also  covered  the 
design  concept  considered  for  transmission  characteristics 
in  the  literature.  Chapter  2 covers  the  force  analysis  for 
the  higher  pairs  including  friction  forces.  A discussion  of 
design  criteria  including  Coulomb  frictional  effects  is 
covered  in  Chapter  3.  Chapter  4 carries  out  the  static 
force  and  torgue  analysis,  which  is  needed  to  define  the 
mechanical  advantages  and  the  instantaneous  efficiency. 

Then,  Chapter  5 gives  the  formulations  of  optimal  design  ac- 
cording  to  the  transmission  criteria  defined  in  the  previous 
chapters  and  other  necessary  constraints.  Three  numerical 
examples  are  given  in  Chapter  6,  which  are  R-Sp-R,  R-Cp-C, 
and  R-Pl-C-R  mechanisms,  respectively.  The  R-Cp-C  mechanism 
is  statically  indeterminate,  and  therefore  it  is  modeled 
with  specified  output  resistant  force  and  torque.  Conclu- 
sions  and  recommendations  for  future  research  are  given  in 
the  final  chapter. 


CHAPTER  2 

THE  MECHANISMS  AND  THE  HIGHER  PAIRS 
2 . 1 The  Mechanisms 

Higher  pairs  can  be  represented  by  combinations  of 
prismatic  and  revolute  pairs  with  special  axis  locations  and 
orientations.  This  means  that  the  direct  constraints  can  be 
formulated  for  the  higher  pairs  themselves.  This  simplifies 
the  synthesis  and  analysis  of  the  mechanisms  with  higher 
pairs.  Due  to  the  multiple  degrees  of  freedom  of  relative 
motion  of  higher  pairs,  the  number  of  links  for  a certain 
function  can  be  reduced  as  compared  to  lower-pair 
mechanisms.  This  means  that  more  compact  mechanisms  can  be 
constructed  by  the  use  of  higher  pairs.  The  kinematic  syn- 
thesis and  analysis  of  mechanisms  containing  higher  pairs 
has  been  shown  by  Ghosal  (1982)  for  three,  four  and  five- 
link  mechanisms.  Displacement  equations  based  on  Ghosal 's 
work  are  listed  in  the  appendices  for  two  of  the  numerical 
examples  in  later  chapters.  The  displacement  equations  for 
one  of  the  examples,  namely  the  R-P^-C-R  mechanism,  are  not 
in  the  existing  literature.  These  equations  are  derived  in 
a similar  manner  as  those  in  Ghosal 's  work  and  are  il- 
lustrated in  Appendix  F.  In  constructing  the  mechanisms,  if 
the  two-DOF  lower  pair  is  grounded  instead  of  being  a float- 
ing joint,  the  mechanisms  so  formed  will  be  statically 
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indeterminate.  For  example,  the  R-Cp-C  and  R-PL-R-C 
mechanisms  are  statically  indeterminate  if  only  the  input 
torque  is  specified. 

2.2  Force  Analysis  for  Higher  and  Lower  Pairs 
Without  Coulomb  friction,  the  force  analysis  is  rather 
straightforward  for  higher  pairs,  especially  for  those 
having  a plane  element,  since  the  reaction  force  is  always 
along  the  direction  of  the  normal  of  the  contact  surface. 
However,  when  Coulomb  friction  is  involved,  the  direction  of 
the  reaction  force  will  be  affected  by  the  direction  of  the 
relative  velocity  between  the  contact  pair  elements.  The 
analysis  procedure  will  be  shown  later. 

2.2.1  Slotted-sphere ( ST| ) pair 

A slotted-sphere  pair  is  a lower  pair  with  two  degrees 
of  freedom  of  relative  rotation  between  pair  elements  and 
may  be  viewed  as  two  perpendicular  revolute  pairs  with  in- 
tersecting axes.  Figure  (2-1)  shows  an  pair  where  the 
inner  sphere  and  the  spherical  shell  are  denoted  as  the  i-th 
and  j-th  link,  respectively.  The  origins  of  the  two  moving 
coordinates  coincide.  The  relative  rotations  between  the 
two  coordinate  systems  are  denoted  as  ©xj_  j and  ©2i  j > respec- 
tively. One  moment  and  three  force  components  are  acting 
between  the  two  pair  elements.  If  the  i-th  link  is  the  in- 
put link  and  the  j-th  link  is  the  output  link,  the  trans- 
mitted moment,  j , must  be  simultaneously  perpendicular  to 
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Figure  2-1  Slotted-sphere ( SL ) pair 
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zi  and  Xj  axes,  i.e.,  in  the  direction  of  y*.  To  count  the 
frictional  effects  induced  by  this  moment,  a force  couple 
consisted  of  and  are  illustrated  at  the  two  assumed 
contact  points  P-j_  and  P2  to  replace  the  moment,  as  shown  in 
the  force  diagram  for  the  SL  pair  in  Figure  (2-2).  The  two 
forces  forming  the  couple  will  cause  two  Coulomb  sliding 
friction  forces,  Ff^  and  Ff2.  The  directions  of  the  two 
friction  forces  are  opposite  to  the  relative  sliding 
velocities  Vsl  and  Vs2,  respectively.  The  directions  of  the 
two  relative  sliding  velocities  must  be  opposite  to  each 
other,  since  the  relative  motion  is  caused  by  the  relative 
rotations  only,  and  will  produce  tangent  velocities  at  the 
two  contact  points.  Thus,  the  directions  of  the  friction 
forces  Ffl  and  Ff2  are  also  opposite  to  each  other.  The 
resultant  friction  moment,  Mf,  must  be  in  the  direction  of 
the  Xj  axis,  i.e.,  normal  to  Yjtj  plane.  The  transmitted 
force,  Fj_j  , can  have  three  components  and  its  line  of  action 
intersects  the  inner  surface  of  the  spherical  shell  at  the 
modelled  contact  point  P.  Also,  there  is  a sliding  friction 
force,  Ff,  caused  by  Fj_j  at  point  P.  The  direction  of  Ff  is 
opposite  to  the  relative  sliding  velocity  which  is  caused  by 
the  two  relative  rotations  at  point  P. 

To  summarize  the  forces  and  moments  at  the  pair,  we  may 
write  the  resultant  force  at  the  center  of  the  sphere,  as 

— ri j = lij  + If  (2-1) 

where 

— s = Vs/|Vs|;  Vs  = 


If  ^ I li  j I Y.s ' 


j x rf  j_  j . 
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Figure  2-2  Force  diagram  of  SL  pair. 
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The  quantity  y.  is  the  friction  coefficient,  a positive 
scalar ; r is  the  radius  of  the  sphere;  ^.ij  is  the  unit  vec— 
tor  of  the  transmitted  force  EVjj  , fjj  = Fjj  / | F± j | ; fi-jj  is 
the  relative  angular  velocity  of  element  i with  respect  to 
element  j . The  resultant  moment  may  be  expressed 

Mrij  = Mi  j + riLij  x Ff  + 

2 [ ( r + t/2  ) z i + qx-jj  x Ffi  (2-2) 

where 

— fl  - ”4 I Mel  I — si ' 

1 1^1 1 = IHc2 I = I Mi j | / ( r + t/2); 

— Sl  = Ysl/|Vsl|; 

—sl  = £li  j x [(r  + t/2)ZjL  + qx-j_  ] . 

The  quantity  t,  a positive  scalar,  is  the  thickness  of  the 
spheric  shell,  which  may  be  neglected;  q is  the  radius  of 
the  inserted  round  bar  on  the  two  sides  of  the  inner  sphere; 
The  direction  of  Mj_j  is  initially  assumed  having  the  posi- 
tive direction  along  y . if  q and  t are  both  neglected. 
Equation  (2-2)  becomes 

Mrij  = Mi  j + rfjj  x Ff  + 2rz±  x Ffl  (2-3) 

Here,  Mrj_j  is  the  resultant  moment  both  the  external  and 
frictional  moments. 

2.2.2  Sphere-slotted-cylinder(Sc)  pair 

Figure  (2-3)  shows  an  Sc  pair  where  the  sphere  and  the 
cylinder  elements  are  connected  to  the  i-th  link  and  the 
3~th  link,  respectively.  Initially,  the  origins  of  the  two 
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Figure  2-3  Sphere-slotted-cylinder ( Sc)  pair 
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moving  coordinate  systems  coincide.  An  Sc  pair  permits  two 
relative  rotations,  0xij  and  ©yi j ; and  a relative  transla- 
tion, lzij ( initially  assumed  to  be  zero),  along  the  direc- 
tion of  the  Zj  axis.  The  Sc  pair  may  be  treated  as  a com- 
bination of  two  revolute  pairs  and  one  prismatic  pair,  where 
the  two  revolute  axes  intersect  on  the  axis  of  the  prismatic 
pair.  One  moment  and  two  force  components  are  transmitted 
between  the  two  pair  elements,  whose  directions  are  shown  in 
Figure  (2-3).  The  transmitted  moment,  Mj_ j , is  in  the  direc- 
tion of  z , which  is  normal  to  the  Xjy^  plane.  If  the  pair 
is  considered  to  be  in  surface-contact,  the  moment  may  be 
replaced  by  a force  couple,  consisting  of  and  Fc2,  ac- 
ting at  the  two  assumed  contact  points,  and  P2.  Figure 
(2-4)  illustrates  the  force  diagram  of  the  Sc  pair.  The 
vector  of  the  force  couple  is  parallel  to  the  Xj  axis  and 
is  passing  through  the  yi  axis.  Point  V1  is  the  intersec- 
tion with  the  normal  surface  of  the  sphere  rod  and  the  edge 
of  the  cylinder  slot.  Point  P2  is  where  Fc2  going  from  the 
center  of  the  inner  sphere  and  passing  through  the  inner 
surface  of  the  cylinder  element.  Thus,  two  friction  forces, 
Ffl  and  Ff2,  are  induced,  whose  directions  are  opposite  to 
the  directions  of  the  relative  sliding  velocities,  Vsl  and 
— s2'  respectively.  Since  these  two  relative  sliding 
velocities  are  parallel  with  the  YjZj  plane,  the  direction 
of  the  vector  of  the  friction  moment,  Mf , must  be  in  a plane 
normal  to  the  YjZj  plane,  i.e.,  collinear  to  Xj  axis.  The 
transmitted  force,  Fjj , has  two  components  in  the  direction 
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Figure  2-4  Force  diagram  of  Sc  pair. 
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of  Xj  and  yj  axes,  respectively.  Thus,  the  transmitted 
force  lies  on  the  Xjyj  plane.  The  contact  point  P is  the 
intersection  with  the  arc  of  contact  of  the  normal  to  the 
axis  of  the  cylinder,  coplanar  with  the  axis  and  the  line  of 
action  of  the  force  acting  on  and  going  through  the  center 
of  the  sphere.  The  forces  and  moments  may  be  summarized  at 
the  center  of  the  sphere.  The  resultant  force  is 

— ri  j = lij  + If  + Ffl  + Ff2  (2-4) 

where 

— f = “V1 1 Hi  j IZs  ' 

Ffl  = -u|Fcl|vsl; 

— f 2 = “^1 Ic2 IZs2 ; 

and 


Zs  = Yo/ |v 


s'  i _s  I ' 


— s = «ij  x rf i:j  + vzij 


Hell  “ I FC2 1 ~ |Mj_j  |/[0.5(r  + h)/cos0xij]; 
Zsl  = Ysl/|Vsl|; 

Ysl  = «ij  X rxj  + vzi j ; 

—s2  = Ys2/ I YS2 | ; 


Ys2  - £li  j x H(r  + h ) /cos0xj_  j ]yj_  - qx  j } + 
ioi]  ' 

The  quantity  r is  the  radius  of  the  sphere;  q is  the  radius 
of  the  sphere  rod;  h is  the  height  of  the  slot  of  the 
cylinder  element;  Vz^j  is  the  relative  linear  velocity  along 
the  cylinder  axis;  and  M^j  is  initially  assumed  in  the 
direction  of  z axis.  The  resultant  moment  may  be  written 
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— ri j = Mij  + rfj_j  x Ff  + rxj  x Ffl  + 

{ [ ( r + h ) /cos©xi j ] Yi  - qxj}  x Ffl  (2-5) 

When  the  force  F*  of  the  sphere  element  exerts  onto  the 
grooved  cylinder  element  is  directed  toward  the  open  slot  in 
the  cylinder,  then  P splits  into  two  points.  Although  the 
exact  contact  points  P-^  and  P2  are  at  the  edges  of  the 
slot  as  shown  in  Figure  (2-5),  we  can  use  the  theoretical 
contact  point  P*,  which  lies  on  the  line  of  action  of  F* 
through  the  center  of  the  sphere,  since  this  does  not  change 
the  resultant  force. 

2.2.3  Sphere-groove ( Sg)  pair 

The  Sg  pair  allows  three  relative  rotations  and  one 
translation  between  the  pair  elements.  An  Sg  pair  may  be 
viewed  as  a combination  of  three  revolute  pairs  and  one 
prismatic  pair,  where  the  three  revolute  pair  axes  are  con- 
current on  the  axis  of  the  prismatic  pair.  Figure  (2-6)  il- 
lustrates an  Sg  pair  where  the  inner  sphere  and  the  groove 
are  denoted  as  the  i-th  and  the  j-th  link,  respectively. 

There  are  two  force  components  transmitted  through  the  two 
pair  elements  and  no  moments  transmitted.  The  location  of 
the  acting  point  P of  the  resultant  force  at  an  Sg  pair  can 
be  described  similarly  to  that  at  an  Sc  pair  (see  Figure 

(2-6)).  The  resultant  force  may  be  summarized  at  the  sphere 
center,  as 

— ri j = lij  + If  (2-6) 

where 
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Figure  2-5  Theoretical  contact  point  of  Sc  or  Sg  pair 
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Figure  2-6  ■=  Sphere-groove (Sg)  pair 
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~VJ- 1 Hi j I ( £li j x rfj_j  + vzij  ) 

Ff ; 

I— i j x rIi  j + Izij  I 

The  quantity  r is  the  radius  of  the  sphere,  is  the  rela- 

tive angular  velocity  vector;  fjj  is  the  unit  vector  of  Fj_j  ; 
— zij  is  the  relative  linear  velocity  along  the  froove.  The 
resultant  moment  at  the  sphere  center  is 

— ri j = rIij  x If  (2-7) 

In  the  force  and  torque  analysis  for  the  mechanisms  contain- 
ing an  Sg  pair,  the  resultant  force  can  be  summarized  at  the 
contact  point  instead  of  the  sphere  center  and  the  resultant 
moment  is  eliminated.  It  is  to  be  noted,  however,  that  when 
the  sphere  rod  reaches  the  edge  of  the  groove,  a moment 
about  the  z*  axis  will  be  transmitted. 

2.2.4  Sphere-plane ( Sp)  pair 

Figure  (2-7)  shows  the  Sp  pair  where  the  inner  sphere 
and  the  plane  elements  are  denoted  as  the  i-th  and  the  j-th 
link,  respectively.  The  pair  permits  three  relative  rota- 
tions and  two  relative  translations.  The  three  rotation 
axes  are  concurrent  at  the  center  of  the  sphere.  There  is 
only  one  force  acting  on  the  contact  surface.  The  contact 
point  P for  a Sp  pair  must  lie  at  the  intersection  of  the 
plane  and  the  force  vector  F^j  passing  through  the  center  of 
the  sphere  and  normal  to  the  plane . The  point  P can  be  on 
either  side  of  the  plane  element  depending  on  the  direction 
of  the  transmitted  force.  When  friction  is  present,  the 
resultant  force  can  be  summarized  at  the  contact  point 
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Figure  2-7 


Sphere-plane ( Sp)  pair. 
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instead  of  the  sphere  center.  Thus,  there  will  be  no  moment 
but  only  forces  acting  at  the  contact  point,  and 

— ri j = lij  + If  (2-8) 

where 

If  “ I li j I Zs 

The  unit  vector  of  the  relative  sliding  velocity  at  the  con- 
tact point,  vs,  can  be  obtained  directly  from  the  displace- 
ment equations. 

2.2.5  Cylinder-planet Cp)  pair 

There  are  two  relative  rotations  and  two  translations 
between  the  pair  elements  of  a Cp  pair  denoted  by  0xij  and 
®zij ' as  shown  in  Figure  (2-8).  In  such  an  interpretation, 
the  axes  of  one  revolute  pair  and  one  prismatic  pair  coin- 
cide, namely  Zj_,  and  the  other  two  pair  axes  are  perpen- 
dicular to  each  other  and  to  the  two  coincident  axes.  The 
relative  translations  are  along  the  yj  and  zj  axes,  respec- 
tively. There  is  only  one  moment  component  and  one  force 
component  transmitted  between  the  pair  elements,  which  are 
acting  on  the  contact  line.  The  contact  line  for  the  Cp 
Pa^r  is  the  intersection  between  the  plane  element  and  the 
normal  plane  that  passes  through  the  axis  of  the  cylinder, 
as  shown  in  Figure  (2-8).  Without  friction,  the  direction 
of  the  transmitted  moment  vector,  , must  be  perpendicular 
to  the  contact  line  and  lies  in  the  contacting  plane  ele- 
ment. The  direction  of  the  transmitted  force,  Fj_j  , must  be 
normal  to  the  contacting  plane  element.  The  moment  vector 
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Figure  2-8  Cylinder-plane ( Cp ) pair 
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Mi j may  be  modeled  as  a force  couple  consisting  of  Fcl  and 
Fc2  acting  at  the  ends  of  the  contact  line  (Figure  (2-9)). 
When  friction  is  present,  two  sliding  friction  forces  are  to 
be  added,  which  lie  on  the  contacting  plane  element  sur- 
faces. Thus  the  direction  of  the  friction  moment,  Mf , must 
be  substantially  normal  to  the  plane  element.  Unlike  the 
previous  pairs,  which  all  have  a sphere  element  and  the 
transmitted  force  components  can  be  regarded  as  a single 
force  passing  through  the  sphere  center,  the  Cp  pair  has  a 
cylinder  element  and  the  transmitted  force  components  are 
distributed  along  the  contact  line  normal  to  the  plane  ele- 
ment. These  forces  will  produce  distributed  friction  forces 
opposing  the  direction  of  the  relative  sliding  velocities, 
which  vary  both  in  direction  and  in  magnitude  along  the  con- 
tact line.  Therefore,  it  is  difficult  to  determine  the 
transmitted  forces  and  moments  analytically  when  considering 
distributed  transmitted  forces. 

To  simplify  the  force  and  torgue  analysis,  the  trans- 
mitted force  and  moment  may  be  modeled  by  two  concentrated 
forces  acting  at  the  two  ends  of  the  contact  line,  where  the 
two  forces  may  have  different  magnitudes  and  their  direc- 
tions are  parallel  to  each  other.  This  is  assuming  that  the 
cylinder  and  the  plane  elements  are  rigid  and  there  may  ex- 
ist some  small  clearance  between  cylinder  and  the  plane  ele- 
ments. As  shown  in  Figure  (2-9),  the  normal  forces  are  Fj_j-j_ 
and  Fij2.  Their  directions  are  initially  assumed  to  be  op- 
posite to  each  other,  parallel  to  the  Xj  axis.  Two  friction 
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Figure  2-9 


Force  diagram  of  Cp  pair. 
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forces  are  induced  at  the  two  contact  points,  and  P2.  To 
determine  the  assumed  location  of  the  resultant  force,  the 
two  normal  forces  are  first  assumed  acting  at  the  axis  of 
the  cylinder  element.  The  forces  may  be  summarized  at  the 
gravity  center  (the  midpoint  of  the  contact  line).  The 
resultant  force  can  be  written 

-rij  = lijl  + Hi j 2 + Ifl  + If 2 (2-10) 

where 

— f 1 = I lijl  llsl/ llsl  I > 

— f 2 = ~U| Fij2 I Vs2/ | Vs2 | ; 

Vsl  = £li  j x 01^  + Vj_  j ; 

Is2  = £li  j x 21 2 + lij  > 

9Rl  = (w/2)z_j_  + qx j ; 

OP2  = ~(w/2)z^  - qx  j . 

The  vector  fij_j  and  V^j  are  the  relative  angular  and  linear 
velocities,  respectively;  OP^  and  OP 2 are  the  vectors  from 
the  midpoint  of  cylinder  axis  in  contact  to  the  two  contact 
points  Px  and  P2,  respectively;  and  q is  the  radius  of  the 
The  transmitted  moment  can  be  written 

-rij  = 2El  x (liji  + Ffi)  + OP2  x ( Fj_  j 2 + Ff  2 ) (2-11) 

If  the  radius  of  the  cylinder,  q,  is  neglected  in  force  and 
torque  analysis,  the  two  vectors  OP]_  and  OP2  can  be 
simplified,  as 

— 1 = ( w/2 ) z.i  (2-12) 

and 

—2  = _(w/2)Zj_,  (2-13) 

respectively.  This  means  that  the  two  normal  forces  and  the 
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two  friction  forces  are  assumed  acting  at  the  axis  of  the 
cylinder. 

2.2.6  Plane (Py)  pair 

A PL  pair  has  planar  surface  contact.  It  is  il- 
lustrated in  Figure  (2-10) . The  pair  allows  one  relative 
rotation,  0xj_ j , and  two  relative  translations,  lyj_j  and 
Izij / between  the  pair  elements.  The  plane  pair  may  be 
viewed  as  a combination  of  one  revolute  pair  and  two  pris- 
matic pairs,  where  the  two  prismatic  pair  axes  intersect  at 
and  the  revolute  pair  axis  is  perpendicular  to  the  two  pris- 
matic axes.  There  are  two  moment  components  and  one  force 
component  acting  between  the  pair  elements.  The  directions 
of  the  transmitted  moments,  Myij  and  Mzij , are  in  the  direc- 
tion of  yj  and  Zj,  respectively;  and  the  direction  of  the 
transmitted  force,  Fj_j  is  normal  to  the  contact  surface, 
i.e.,  in  the  direction  of  Xj_  and  Xj  . Similarly  to  the  Cp 
pair,  the  transmitted  force  is  distributed  on  the  contact 
surface  and  can  not  readily  be  represented  by  a concentrated 
force  and  moment. 

Although  the  contact  in  the  P^  pair  is  a surface  con- 
tact, we  may  assume  three  concentrated  normal  force,  Fiji, 

— i j 2 , and  Iij3'  applied  at  three  of  the  four  corners  of  the 
contact  plane,  as  shown  in  Figure  (2-11).  This  is  also 
based  on  the  assumption  of  the  existence  of  clearance  be- 
tween pair  elements.  It  is  reasonable  to  neglect  the  thick- 
ness of  the  inner  plane  element,  and  the  three  transmitted 
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Figure  2-10 


Plane  (PL)  pair. 


39 


Figure  2-11 


Force  diagram  of  PL  pair. 
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forces  can  be  assumed  acting  at  the  same  contact  surface. 
The  directions  of  the  three  transmitted  forces  are  in  the 
direction  of  the  Xj  axis.  The  three  normal  forces  will 
produce  three  sliding  friction  forces  lying  in  the  contact 
plane.  The  resultant  force  and  moment  may  be  summarized  at 
the  gravity  center  as 
3 

-rij  = jJ-^-ijk  + — fk)  (2-14) 

and 

3 

-rij  = k=x— ^ X (-ijk  + Ifk)  (2-15) 


CHAPTER  3 

TRANSMISSION  CRITERIA 

The  transmission  criteria  for  spatial  mechanisms  can  be 
fully  evaluated,  by  the  deviation  angles  of  each  joint  for 
force  and  moment  transmission,  the  mechanical  advantage  and 
the  instantaneous  efficiency  of  the  mechanism,  and  the  ac- 
celerations for  better  dynamics  performance.  The  following 
sections  define  these  criteria  according  to  the  geometry  of 
each  pair  mentioned  in  Chapter  2 and  for  some  typical 
mechanisms . 


3 . 1 Deviation  Angles 

To  explore  the  transmission  characteristics  of 
mechanisms,  not  only  must  the  force  and  moment  transmission 
from  the  coupler  to  the  output  link  to  be  evaluated,  but 
also  the  force  and  moment  transmission  in  each  floating 
joint.  Good  force  and  moment  transmission  in  each  floating 
joint  provides  the  basis  of  good  transmission  characteris- 
tics for  the  mechanism.  Furthermore,  it  ensures  the 
avoidance  of  separation  of  force-closed  pair  elements.  It 
iJ  therefore  necessary  to  separately  define  deviation  angles 
for  the  internal  force  and  moment  transmission  of  each 
floating  joint  and  for  the  force  and  moment  transmission  at 
the  joint  connecting  the  coupler  link  and  the  output  link 
for  the  output  motion.  In  order  to  define  the  deviation 
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angles,  the  transmitted  force  and  moments  for  the  Cp1  and 

Pairs  regarded  as  acting  at  the  center  of  gravity  of 
the  pair  contact.  When  considering  the  deviation  angles  for 
force  and  moment  transmission,  the  directions  of  the  trans- 
mitted force  and  moment  vectors  are  of  concern. 


3,1,1  Deviation  angles  for  producing  output  force  and  moment 
Figure  (3-1)  illustrates  a coupler  link  and  an  output 
link,  which  are  assumed  to  be  connected  by  a joint,  Q,  that 
wiii  transmit  both  force  and  moment.  This  floating  joint 
may  be  either  a lower  or  higher  pair.  The  grounded  output 
joint  is  assumed  to  be  a cylindrical  joint  so  there  are  both 
angular  and  linear  output  displacements.  Then,  the  devia- 
tion angle  for  the  force  transmitted  to  produce  the  transla- 
tional motion  of  the  output  link  can  be  defined,  as 


6 ft  = cos 


-1 


—CO  * Yot 


(3-1) 


I Yco I I Yot I 

where  Fco  is  the  force  exerted  from  the  coupler  link  onto 
the  output  link  through  the  floating  joint  Q,  Vot  is  the  ab- 
solute velocity  at  the  joint  Q caused  by  the  translational 
motion  of  the  output  joint.  The  deviation  angle  for  the 
force  transmitted  to  produce  the  rotational  motion  of  the 
output  link  is  defined 
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(3-2) 
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Figure  3 1 Deviation  angles  for  producing  output  force  and 
moment . 
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where  Vor  is  the  absolute  velocity  at  the  joint  Q caused  by 
the  rotational  motion  of  the  output  joint.  The  deviation 
angle  for  the  moment  transmitted  to  produce  the  rotational 
motion  of  the  output  joint  is 

_i  —co  ‘ Ho 

6MR  “ cos  (3-3) 

I Moo  I 

where  is  the  moment  exerted  by  the  coupler  to  the  output 

link  at  the  joint  Q;  and  Uq  is  the  unit  vector  of  the  direc- 
tion of  the  output  joint  axis,  i.e.,  the  direction  of  the 
output  resistant  moment  vector.  For  the  four-link 
mechanisms,  the  moment  exerted  by  the  coupler  link  to  the 
output  link  at  the  Q joint  may  consist  of  two  components, 
one  is  caused  by  the  driving  torque  and  the  other  one  is 
caused  by  the  transmitted  force  at  the  previous  floating 
joint.  Thus,  for  the  four-link  mechanisms,  there  will  be 
two  deviation  angles  to  distinguish  the  effort  of  each  mo- 
ment. Similarly,  for  five-link  mechanisms  the  exerted  mo- 
ment at  the  Q joint  may  consist  of  three  components,  one  of 
them  is  caused  by  the  driving  torque  directly  and  the  other 
two  components  are  caused  by  the  two  transmitted  forces  at 
the  previous  two  floating  joints,  respectively.  If  the 
cosine  of  these  deviation  angles  approaches  to  zero,  there 
may  be  redundant  freedoms  generated  in  the  output  link, 
i.e.,  the  transmitted  force  or  moment  does  not  contribute  to 
the  output  force  or  moment.  For  example:  cos6FT  = 0 indi- 
cates there  is  no  translation  displacements  at  the  output 
joint;  both  cos6fr  = o and  cos6MR  = 0 indicate  that  the 
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rotational  displacements  is  lost  in  the  output  joint.  In 
practice,  the  absolute  value  of  the  cosine  of  the  minimum 
deviation  angle  may  be  used  to  define  the  design  criteria, 
as 

| cosSj^h  | > (3  (3-4) 

where  (3  is  a positive  value  less  than  one  and  away  from 
zero . 

Eguations  (3-1)  to  (3-3)  apply  to  any  pair  that  is  ad- 
joining the  coupler  link  and  the  output  link.  These  three 
deviation  angles  are  for  the  use  of  defining  how  the  exerted 
force  and  moment  from  the  coupler  link  to  the  output  link 
affect  the  resistant  force  and  moment  at  the  output  joint. 


^ ^ • 2 Deviation  angles  for  internal  force  and  moment 
transmission ——  ~ 

In  cam  mechanisms,  the  force  transmission  characteris- 
tic between  the  cam  and  the  follower  is  defined  by  the  pres- 
sure angle,  which  is  the  angle  between  the  direction  of  mo- 
tion of  the  follower  and  the  normal  to  the  cam  surface  at 
the  contact  between  the  cam  and  the  follower  (Chen,  1982). 

In  higher  pairs  other  than  the  cam  pair,  a similar  defini- 
tion can  be  applied. 

For  either  the  higher  or  the  lower  pair  that  is  located 
in  the  joint  connecting  the  coupler  link  and  the  output 
link,  the  angles  between  the  directions  of  the  transmitted 
force  and  moment  vectors  at  the  contact  point  or  center  of 
gravity  of  the  contact  line,  curve,  or  area  and 
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the  direction-of -motion  vectors  serve  as  the  deviation 
angles.  The  friction  angle  is  also  involved  in  determining 
the  transmission  criteria.  The  deviation  angles  measure  how 
the  coupler  force  and  moment  produce  the  linear  and  angular 
motion  at  the  output  joint.  This  effect  upon  the  transmis- 
sion characteristics  between  the  coupler  and  the  output  link 
has  been  shown  in  the  previous  section.  Besides  these 
deviation  angles,  there  are  also  deviation  angles  involving 
internal  force  and  moment  transmission  needs  to  be  con- 
sidered between  other  mechanism  joints  especially  when  the 
mechanism  contains  a higher  pair.  This  is  due  to  the  nature 
of  the  external-groove,  cylinder  and  plane  elements  of 
higher  pairs  in  transmitting  force  and/or  moment.  When  the 
forces  are  transmitted  through  an  Sg3  or  Sc4  pair,  the 
transmitted  forces  must  lie  on  the  plane  that  contains  the 
contact  point  and  perpendicular  to  the  axis  of  the  external 
grove  or  the  cylinder  element  as  illustrated  in  Chapter  2. 
Thus,  the  direction  of  the  transmitted  force  is  restricted, 
as  is  also  the  transmitted  moment.  If  the  deviation  angle 
for  the  force  or  moment  transmission  approaches  to  90  de- 
grees, the  higher  pairs  may  be  separated  or  generate  rela- 
tive passive  freedom  in  the  higher  pair.  Similarly,  for  the 
mechanisms  containing  a Sp5,  Cp,  or  PL  pair  which  have  an 

3.  Sphere-in-groove 

4.  Sphere-in-cylinder 

5.  Sphere-on-plane 
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external  plane  element,  the  transmitted  force  must  be  along 
the  direction  of  the  normal  of  the  contact  plane  and  the 
direction  of  the  vector  of  the  transmitted  moment  must  be  on 
the  contact  plane.  Even  in  lower  pairs  such  as  the  revolute 
(R),  prismatic  (P),  or  the  cylindrical  (C)  joint,  the  direc- 
tion of  the  transmitted  force  and  moment  may  also  be 
restricted.  Therefore,  it  is  necessary  to  discuss  the 
definition  of  deviation  angles  for  internal  force  and  moment 
transmission  according  to  the  geometrical  constraints  of 
each  pair. 

3. 1.2.1  SL( slotted  sphere)  pair 

The  direction  of  the  force  transmitted  between  the  pair 
elements  of  an  SL  pair  is  not  restricted  since  it  does  not 
have  any  relative  translational  freedoms.  When  an  SL  pair 
is  connecting  the  coupler  link  and  the  output  link,  Equation 
(3-1)  to  (3-2)  for  force  transmission  can  be  applied. 
However,  the  direction  of  the  transmitted  moment  is 
restricted,  which  has  been  shown  in  Chapter  2.  Thus,  a 
deviation  angle  for  the  internal  moment  transmission  needs 
to  be  defined.  An  R-SL  dyad  is  shown  in  Figure  (3-2),  where 
the  R joint  is  assumed  to  be  the  input  joint  with  an  input 
driving  torque  M^.  The  inner  sphere  is  connected  to  the  in- 
put link.  The  input  moment  is  applied  from  the  inner  sphere 
to  the  outside  sphere;  the  direction  of  the  transmitted  nor- 
mal moment  vector  M23  must  be  along  the  direction  of  the  y* 
axis,  which  lies  in  the  y3z3  plane,  i.e.,  the  plane  in 
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Figure  3-2  Deviation  angles  of  an  SL  pair. 


^7/ 
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which  the  centerline  of  the  cut  slot  of  the  external  sphere 
lies.  The  deviation  angle  for  the  moment  transmission  may 
therefore  be  expressed: 

6M23  = cos_1u2  • y*  (3-5) 

where  u2  is  the  unit  vector  of  the  input  joint  axis  and  y* 
is  the  unit  vector  resulted  by  rotating  y2  about  the  x2  axis 
by  the  angle  of  0x22. 

3. 1.2. 2 Cp  and  Sp  pairs 

For  a three- link  mechanism  containing  an  Sp  or  Cp  pair, 
the  only  floating  joint  is  the  higher  pair,  where  the  trans- 
mitted forces  must  be  along  the  direction  of  the  normal  of 
the  contacting  tangential  plane  and  the  transmitted  moment 
must  lie  on  the  contacting  plane.  Thus,  the  deviation  an- 
gels for  internal  force  and  moment  transmissions  are  ex- 
plicitly defined  by  the  transmitted  normal  force  and  moment 
and  the  directions  of  the  output  velocity  which  causes  the 
translation  or  rotation  of  the  output  joint.  However,  this 
is  assuming  that  the  plane  element  is  connected  to  the  input 
link.  If  the  plane  element  is  attached  to  the  output  link, 
then  the  deviation  angles  are  defined  otherwise.  Figures 
(3-3)  and  (3-4)  illustrate  the  two  different  constructions 
of  R-Cp  dyads.  In  Figure  (3-3),  the  cylinder  element  is 
connected  to  the  input  link  denoted  as  link  2.  Two  moving 
coordinates  are  attached  at  the  intersection  of  the  two  con- 
tact forces.  The  deviation  angle  for  force  transmission, 
®F23'  is  the  angle  between  the  unit  vector  of  the  direction 
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Figure  3-3  Deviation  angles  of  a Cp  pair  (the  cylinder 
element  is  connected  to  the  input  link). 
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Figure  3-4  Deviation  angles  of  a Cp  pair  (the  plane  element 
is  connected  to  the  input  link). 
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of  the  input  effort  force,  F2,  caused  by  the  input  torque  at 
the  center  of  gravity  of  the  line  contact  and  the  direction 
of  the  transmitted  force  normal  to  the  plane  element,  £23. 
Thus,  we  have 


where  Vp  is  the  unit  vector  of  the  absolute  velocity  Vp  and 
X3  is  the  unit  normal  of  the  contact  plane  in  the  direction 
from  the  cylinder  element  to  the  plane  element.  If  link  3 
is  the  output  link,  such  as  in  an  R-Cp-C  mechanism,  the 
deviation  angle  between  the  transmitted  normal  force  and  the 
absolute  translational  motion  caused  by  link  3 is  predeter- 
mined by  the  geometry  of  link  3,  i.e.,  the  value  determined 
by  Equations  (3-1)  is  fixed.  The  other  variable  deviation 
angle  for  force  transmission  is  the  angle  between  the  normal 
force  vector  and  the  absolute  rotational  output  motion.  It 
can  be  defined  as 


where  Vpi  is  the  unit  vector  of  the  absolute  velocity  Vpi 
starting  out  from  the  output  joint.  The  deviation  angle  for 
internal  moment  transmission  is  defined 


where  u2  is  the  unit  vector  of  the  input  joint  axis  and  y* 
is  the  unit  vector  of  the  direction  of  the  moment  vector 
transmitted  from  the  cylinder  element  to  the  plane  element. 
If  link  3 is  the  output  link,  the  deviation  angle  between 
the  transmitted  moment  vector  and  the  absolute  rotational 
motion ( angular  velocity  vector)  is  also  a variable 


6f23  = cos  Xvp  • x3 


(3-6) 


(3-7) 


6M23  = cos  1li2  • 1 


(3-8) 
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parameter,  which  may  be  defined  by  using  Equation  (3-3). 

On  the  other  hand,  when  the  plane  is  connected  to  the 
input  link,  as  shown  in  Figure  (3-4),  the  deviation  angles 
6F23  and  6fr  are  defined  similarly  to  Equations  (3-6)  and 
(3-7).  But,  the  deviation  angle  for  producing  translational 
output  motion  is  variable  now.  It  can  be  defined  as 

6ft  = cos-1X2  * H3  (3-9) 

where  X2  is  the  unit  vector  of  the  normal  of  the  contact 
plane  and  u.3  is  the  unit  vector  of  the  output  joint  axis. 

The  deviation  angles  for  moment  transmission  are  the  same  as 
the  in  previous  case. 

It  can  be  observed  from  Figure  (3-3)  when  6F23  ap- 
proaches to  90  degrees,  the  relative  sliding  between  the 
cylinder  and  the  plane  will  approach  a maximum.  This  is  the 
so  called  indeterminate  slipping,  which  will  cause  the  un- 
certainty configuration  as  defined  by  Hunt  (1979).  Also, 
when  Sjyj2 3 goes  to  90  degrees,  there  will  be  a rotational 
redundant  freedom  between  the  two  pair  elements.  Similar 
definition  applies  to  the  four,  or  five-link  mechanisms  that 
have  a higher  pair  with  a plane  element. 

For  the  Sp  pair,  as  shown  in  Figure  (3-5),  the  defini- 
tion for  the  internal  force  transmission  is  similar  to  that 
for  the  Cp  pair,  except  that  the  Sp  pair  does  not  transmit 


moment . 
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Figure  3-5  Deviation  angle  of  an  Sp  pair. 
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3 . 1 . 2 . 3 PL  pair 

For  the  PL  pair,  the  deviation  angle  for  the  internal 
force  transmission  is  a variable  parameter  since  the  posi- 
tions where  the  forces  applied  are  variable.  Figure  (3-6) 
shows  an  R-PL  dyad,  where  the  R joint  is  an  input  joint  with 
an  input  moment  M^.  The  deviation  angle  6F23  is  the  angle 
between  the  input  effort  force  F2  anc^  the  transmitted  force 
£23*  The  transmitted  moment  vector  must  lie  on  the  contact 
surface.  As  shown  in  Chapter  2,  the  moment  has  two  com- 
ponents on  the  Yj  and  Zj  axes,  respectively.  The  direction 
of  the  transmitted  moment  vector  M23  is  the  projection  of 
the  driving  moment  vector  on  the  contact  surface.  Thus, 
the  deviation  angle  6^23  for  the  moment  transmission  is 
defined  as  the  angle  between  the  vectors  of  the  driving 
torque  and  the  transmitted  moment.  6^2 3 is  a fixed  value 
here . 

3. 1.2. 4 Sg  and  Sc  pairs 

The  deviation  angles  of  force  transmission  for  the  Sg 
and  Sc  pairs  are  defined  similarly  as  those  for  the  Sp  and 
Cp  pairs,  except  that  the  direction  of  the  transmitted  force 
is  now  located  on  the  plane  perpendicular  to  the  joint  axis 
and  passes  through  the  center  of  the  sphere  and  the  contact 
point.  Figure  (3-7)  illustrates  an  R-Sc  dyad  where  the 
sphere  is  connected  to  the  link  2 with  an  R input  joint  and 
the  cylinder  element  is  connected  to  link  3.  Two  coor- 
dinates are  located  at  the  sphere  center  as  described  in 
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Figure  3-6  Deviation  angles  of  a PL  pair 
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Figure  3-7  Deviation  angles  of  a Sc  pair. 
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Chapter  2.  The  deviation  angle  for  the  force  transmission 
is  defined,  as 

6f23  = cos_1vp  • x*  (3-10) 

where  vp  is  the  direction  of  the  input  effort  force,  i.e., 
the  unit  vector  of  the  absolute  velocity  of  the  contact 
point;  and  x~  is  the  unit  vector  of  the  transmitted  force, 
which  is  the  projection  of  vp  on  the  plane  passing  through 
the  sphere  center  and  normal  to  the  axis  of  the  cylinder. 

The  deviation  angle  for  the  moment  transmission  is 

6M23  = cos-1u2  • z*  (3-11) 

where  u2  is  the  direction  of  the  unit  vector  of  the  driving 
torque,  and  z*  is  the  direction  of  the  transmitted  moment 
vector,  which  lies  on  the  y3  and  z3  plane.  Similar  defini- 
tion applies  to  the  mechanisms  containing  an  Sg  pair  except 
there  is  no  moment  transmission. 

3. 1.2. 5 R,  P,  and  C pairs 

The  revolute  and  prismatic  pairs  do  not  have  moment  and 
force  transmitted  along  the  direction  of  the  joint  axis, 
respectively.  There  is  neither  moment  nor  force  transmitted 
along  the  joint  axis  of  a cylindric  pair.  These  restraints 
of  the  force  and  moment  for  the  lower  pairs  also  affect  the 
transmission  characteristics  of  the  mechanisms  since  the 
force  or  moment  component  in  the  direction  of  the  joint  axis 
vanishes.  Thus,  the  deviation  angles  for  the  internal  force 
and  moment  transmission  of  lower  pairs  should  also  be  of 
concern.  To  define  the  deviation  angels  for  these  three 
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lower  pairs,  an  R-C  dyad  is  shown  in  Figure  (3-8),  where  a 
force  and  a moment  are  exerted  at  the  C joint,  which  are  in- 
duced by  the  driving  torque  at  the  R joint.  The  joint  axes 
of  the  R and  C joint  are  denoted  as  u.j_  and  Uj  , respectively. 
The  transmitted  force  and  moment  may  be  expressed 

-ij  = Hj  x £i}  x Uj  (3-12) 

and 


-ij  “ ^ — j x ?2d)  x Uj  (3-13) 

respectively.  Fj_  is  the  effort  force  caused  by  the  input 
torque  at  the  R joint  applied  to  the  j-th  link  through  the  C 
joint,  which  is  in  the  direction  of  the  absolute  velocity  of 
the  i-th  link  at  the  C joint,  v-^;  and  M^  is  the  driving 
torque  vector,  which  is  in  the  direction  of  the  joint  axis, 

— i * Then,  the  cosine  value  of  the  deviation  angles  for  in- 
ternal force  and  moment  transmission  of  the  C joint  can  be 
expressed 


cos6 


Fij 


£i  ' lij 

lEillFii 


(3-14) 


and 


cos6 


'Mi  j 


^i 


— i j 


(3-15 


|M±| iMij I 

The  two  values  of  Equations  (3-14)  and  (3-15)  are  equal  to 
the  sinoFj_ j and  sinoMj_j  , respectively,  where 


7Fi  j 


= cos 


-i 


±1 


Hi  I 


(3-16) 


and 
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Figure  3-8  Deviation  angles  of  a lower 


pair. 
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-1  -1  -3 

°Mij  - cos  — ; ; (3-17) 

I Mi  | 

Equations  (3-16)  and  (3-17)  also  apply  to  the  P or  R joint 
in  place  of  the  C joint. 

3 • 2 Force  and  Moment  Transmissivities 
1^  the  previous  sections,  the  deviation  anqles  for 
producing  output  force  and  moment  and  for  the  internal  force 
and  moment  transmission  of  both  higher  pairs  and  lower  pairs 
were  introduced.  Using  these  deviation  angles  in  practice, 
the  force  and  moment  transmissivities  are  now  to  be  defined. 
This  will  expand  the  usefulness  of  the  deviation  angles. 

The  force  and  moment  transmissivities  are  defined  as  the  ef- 
fectiveness of  how  the  input  effort  force  and  input  moment 
affect  the  output  resistance  force  and  moment.  There  are 
four  force  and  moment  transmissivities:  force-to-f orce , 
force-to-torque,  torque-to-torque , and  force  torque-to- 
torque.  The  force-to-f orce  transmissivity  is  defined  as  how 
the  driving  force  of  an  input  prismatic  joint  or  input  ef- 
fort force  caused  by  an  input  revolute  joint  at  the  joint 
connecting  the  input  link  and  the  coupler  link  affects  the 
translational  motion  of  the  output  link,  or  in  other  words, 
how  it  produces  the  output  force.  The  force-to-torque 
transmissivity  is  defined  as  how  the  driving  force  or  the 
input  effort  force  affects  the  rotational  motion  of  the  out- 
put link  or  how  it  produces  the  output  torque.  The  torque- 
to-torque  transmissivity  is  how  the  driving  torque  intro- 
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duced  at  an  input  revolute  joint  or  the  input  effort  torque 
caused  by  the  driving  force  at  an  input  prismatic  joint  or 
at  the  joint  connecting  the  input  link  and  the  coupler  link 
affects  the  rotational  motion  of  the  output  link  or  produces 
the  output  torque.  For  f our-or-more-link  mechanisms,  the 
exerted  force  at  a floating  link  will  produce  a moment  at  a 
successive  floating  joint.  The  f orce-torque-to-torque 
transmissivity  is  defined  as  how  this  moment  produces  the 
output  torque.  These  transmissivities  are  determined  by  the 
geometry  of  the  mechanism.  The  following  sections  will  dis- 
cuss the  usage  of  these  transmissivities  according  to  the 
construction  of  the  mechanisms. 

3.2.1  R-Sp-R  mechanism 

Figure  (3-9)  shows  an  R-Sp-R  mechanism  where  the  inner 
sphere  is  connected  to  the  input  link.  The  contact  point  is 
denoted  as  P;  P and  P'  are  the  position  vectors  from  the  in- 
put joint  axis  and  the  output  joint  axis  to  the  contact 
point,  respectively;  a moving  coordinate  system  is  attached 
to  the  contact  point,  which  consists  of  the  orthonormal  unit 
vectors  a,  b,  and  c,  where  a is  the  normal  of  the  plane  ele- 
ment. The  force-to-torque  transmissivity  for  mechanisms 
having  such  geometry  can  be  defined,  as 

f FM  = | cos6F23COs6FR|  (3-18) 

where  6p23  defined  similarly  to  Equation  (3-6),  i.e.,  the 
deviation  angle  between  the  input  effort  force,  £2,  and  the 
transmitted  normal  force,  F23/"  and  <5FR  is  defined  similarly 
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Figure  3 9 Transmissivities  of  an  R-Sp-R  mechanism. 
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to  Equation  (3-2),  i.e,  the  angle  between  the  direction  of 
the  transmitted  normal  force,  F23,  and  the  direction  of  the 
output  joint  resistant  force (or  the  absolute  velocity  of  the 
contact  point  as  a point  of  the  output  link).  The  deviation 
angles  of  Equation  (3-18)  can  be  expressed  in  terms  of  the 
geometrical  parameters,  as 

c°s6F23  = vp  • a (3-19) 

and 

cos6fr  = a • yp.  (3-20) 

where  Vp  and  Vpi  are  the  unit  vectors  of  the  absolute 
velocities  of  the  contace  point  as  a point  of  the  input  or 
output  link,  respectively;  a is  the  unit  vector  of  the  nor- 
mal of  the  plane  element. 

3.2.2  R-Cp-C  mechanism 

For  the  mechanisms  having  geometry  as  the  R-Cp-C 
mechanisms  shown  in  Figure  (3-10) , there  are  two  force 
transmissivities,  since  the  output  joint  has  both  transla- 
tional and  rotational  motions.  There  is  also  one  moment 
transmissivity.  The  center  of  gravity  of  the  contact  line 
is  denoted  as  P and  the  normal  of  the  plane  element  is 
denoted  as  a;  b is  the  unit  vector  on  the  contact  plane 
perpendicular  to  a;  the  axis  of  the  cylinder  element  is  c*. 
The  force- to- torque  transmissivity  to  produce  rotational  mo- 
tion of  the  output  link  at  the  output  joint  is  defined 
similarly  to  Equation  (3-18).  The  additional  force-to-force 
transmissivity  is  defined  as 
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Figure  3-10  Transmissivities  of  an  R-Cp-C  mechanism. 
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fFF  = |cos6f23cos6ft|  (3-21) 

where  6F23  is  defined  similarly  to  Equation  (3-19);  and  6FT 

is  defined  similarly  to  Equation  (3-1)  and  can  be  determined 
by: 

cos6ft  = a • u3  (3-22) 

where  u3  is  the  unit  vector  of  the  output  joint  axis.  Also, 
there  is  a torque-to-torque  transmissivity  for  the  input  mo- 
ment to  produce  the  rotational  motion  of  the  output  link  at 
the  output  joint,  which  can  be  expressed 

^MM  - | cos6]y[23COs6[y[R  | (3-23) 

where  6M23  is  defined  by  Equation  (3-8),  which  is  the  angle 
between  the  driving  torque  vector  and  the  transmitted  normal 
moment  vector;  6MR  is  defined  similarly  to  Equation  (3-3), 
which  is  the  angle  between  the  transmitted  normal  moment 
vector  and  the  output  resistant  moment  vector.  The  cosine 
of  deviation  angles  of  Equation  (3-23)  can  be  expressed  by 
the  geometrical  parameters  disregarding  the  sign,  as 

c°s6M23  = H2  ’ k*  (3-24) 

and 

cosSiyiR  = b*  • u3  (3-25) 

where  u2  and  u3  are  the  unit  vectors  of  the  input  and  output 
joint  axes,  respectively.  Unlike  the  deviation  angles  for 
the  force  transmission,  the  values  of  cos6M23  and  cos6MR  are 
variable  parameters,  since  the  direction  of  b*  is  determined 
both  the  axis  of  the  cylinder  element  and  the  normal  of  the 
plane  element,  i.e.,  b*  = a x c*. 
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3.2.3  R-Py -C-R  mechanism 

For  four  or  more-link  mechanisms,  the  deviation  angles 
at  all  floating  joints  are  of  concern.  Consider  an  R-PL-c-R 
mechanism  as  shown  in  Figure  (3-11).  An  input  effort  force, 
Fi,  and  moment,  M^,  are  caused  by  the  input  driving  torque 
at  the  PL  joint.  We  may  assume  the  transmitted  normal  force 
£23  is  in  the  same  direction  as  the  normal  of  the  plane 
elemnt  of  the  P-^  pair,  a,  since  the  sign  of  the  cosine 
values  are  not  concerned.  The  force- to-torque  transmis- 
sivity is  defined 

fFR  = !cos6F23COs6F34cos6FR|  (3-26) 

where  6F23  is  the  angle  between  the  input  effort  force  and 
the  normal  of  the  plane  elements;  6F34  is  the  angle  between 
the  coupler  force  and  the  transmitted  force,  F34>  at  the 
floating  C joint.  The  transmitted  force  F34  can  be  ex- 
pressed, as 

£34  = <U3  x £23 > x u3  (3-27) 

5F34  is  the  angle  between  the  normal  of  the  plane  element 
and  the  unit  vector  of  the  direction  of  F34,  whose  consine 
value  is  the  same  as  the  sine  value  of  the  angle  contained 
by  £23  and  u3.  Thus,  we  may  write 

cos6F34  = sin(cos~1a  • u3 ) (3-28) 

It  can  be  observed  that  6F34  is  a fixed  angle  since  the 
plane  element  is  attached  to  joint  3.  The  third  deviation 
angle  6FR  is  defined  similarly  to  Equation  (3-2),  i.e.,  the 
angle  between  £34  and  the  direction  of  the  resistant  force 
for  the  rotation  of  the  output  link,  which  can  be  expressed 
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Figure  3-11  Transmissivities  of  an  R-PL-C-R  mechanism. 
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cos6F£  — ^34  • V3 1 

(u.3  x a)  x u3  • V31 

= ; (3-29) 

I ( H3  x a)  x u3  | 

where  v3 , is  the  unit  vector  of  the  direction  of  the  ab- 
solute velocity  of  the  center  of  the  floating  C joint  as  a 
point  of  the  output  link.  The  torgue-to-torque  transmis- 
sivity is 

fMM  = !cos%23cos%M34cos6MrI  (3-30) 

The  fixed  angle  6^23  is  the  deviation  angle  of  the  internal 
moment  transmission  of  the  joint.  It  is  defined  by 

cos6M23  = u 2 * 0123  (3-31) 

where 

(a  x u2)  x a 

m23  = — 

I ( a x u2 ) x a | 

is  the  unit  vector  of  the  direction  of  the  transmitted  mo- 
ment vector  M23«  Equation  (3-31)  can  be  also  be  written,  as 
c°s6M23  = sin(cos-1u2  • a)  (3-32) 

The  angle  6^34  is  the  deviation  angle  for  the  sequential 
moment  transmission  at  the  floating  R joint,  i.e.,  the  angle 
between  the  vectors  of  M23  and  M3 4.  The  unit  vector  of  the 
direction  of  the  transmitted  moment  through  joint  3 can  be 
expressed  as 

(u3  x m22 ) x u.3 

~34  = ~777 I T (3-33) 

I (1±3  x S23)  x )±3  I 

Then,  cos6jv[34  can  be  written 
cos6jyj34  = m23  • m3  4 

= sin( cos-1m23  * a3 ) 


(3-34) 
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The  third  angle  6MR  is  defined  as  the  angle  between  the 
transmitted  moment  vector  at  the  floating  C joint  and  the 
output  resistant  torque  vector,  i.e., 

(u3  x m23)  x u3  • u4 

cos6MR  = ; (3-35) 

I (u.3  x m23 ) x u3| 

The  transmitted  normal  force  at  the  PL  joint  will 
produce  a moment  about  the  floating  C joint.  Then,  the 
force  torque-to-torque  transmissivity  for  this  moment  to 
produce  the  output  torque  is 

fpMM  = lc°s6F23COs6FM34cos6FMR|  (3-36) 

The  angle  of  6F23  is  defined  similarly  to  Equation  (3-19). 
The  direction  of  the  force-produced  moment  can  be  written 

— F23  = £3  x a / |P3  x a | (3-37) 

where  P3  is  the  distance  vector  from  joint  3 to  the  center 
of  gravity  of  the  contact  surface  of  the  PL  pair.  The  unit 
vector  of  the  direction  of  the  transmitted  moment  at  the  C 
joint  is 

— FM3 4 = (u3  x mF23)  x u3  / | (u3  x mF23)  x u3 1 
[u.3  x (P3  x a)  ] x u.3 

— — — — ( 3_3g  ) 

I £3  x a| | [u3  x (£3  x a) ] x u3 | 

Thus,  cos6F]v[34  can  written,  as 

cos6fm34  = mF23  • mFM34 

= sin(cos“l(u.3  • P3  x a)  (3-39) 

The  angle  6F]yiR  is  the  angle  between  the  transmitted  moment 
and  the  output  resistance  torque. 

cos6FMR  = Q)FM3  4 * H4 


It  can  be  defined  as 
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[u3  x ( P3  x a) ] x u3  • u4 

TTT I (3-40) 

I —3  x aj  | [u3  x (P3  x a) ] x u3 | 

3 . 3 Frictional  Effects 

In  Figure  (3-12),  a sphere  and  a plane  are  in  contact 
at  point  P,  where  the  sphere  is  connected  to  the  input  link 
which  has  revolute  joint  and  the  plane  is  connected  to  the 
output  link  which  has  the  other  revolute  joint.  The  ab- 
solute velocity  at  the  contact  point  P as  a point  of  the  in- 
put link  is  Vp  and  as  a point  of  the  output  link,  it  is  V 1 . 
Then,  the  relative  velocity  is  denoted  by  vnn,  where  V , = 
-p  ~ Ypi . The  sliding  relative  velocity,  i.e.,  the 
projected  component  of  the  relative  velocity  on  the  contact 
plane  is  given  as 

-s  - -PP'~  (-PP*  * (3-41) 

where  a is  the  unit  vector  of  the  normal  of  the  contact 

plane.  Now,  if  there  is  an  input  torque  applied  at  the  input 
link  at  the  input  revolute,  an  input  effort  force  f2  , which 
is  in  the  direction  of  Vp,  acts  on  the  contact  point  from 
the  sphere  to  the  plane.  The  input  effort  force  will  cause 
a normal  force  Fn  applied  from  the  sphere  to  the  plane 
through  the  contact  point.  The  deviation  angle  in  this  case 
then  is  defined  as  the  angle  between  F2  and  Fn.  If  the 
roles  of  the  sphere  and  the  plane  are  reversed,  the  devia- 
tion angle  needs  to  be  defined  differently,  as  discussed  in 
the  previous  sections.  Also,  there  is  another  force  com- 
ponent of  the  input  effort  force  projected  on  the  contact 
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Figure  3 12  Frictional  effects  in  plane  motion. 
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plane,  which  is  denoted  as  Fs.  This  is  the  force  causing 
the  relative  movement  between  the  sphere  and  the  plane  on 
the  contact  plane.  Now,  if  there  is  friction  between  the 
sphere  and  the  plane,  a friction  force,  Ff,  will  be  induced. 
The  direction  of  the  friction  force  is  opposite  to  the 
direction  of  the  relative  sliding  velocity,  Vs,  as  shown  in 
Figure  (3-12).  The  resultant  reaction  force,  Fr,  consists 
of  the  resistant  normal  force  and  the  friction  force.  It 
can  be  observed  that  the  friction  force  must  be  overcome  by 
the  sliding  force  in  order  to  keep  the  mechanism  in  motion. 
This  can  be  easily  followed  by  first  considering  the  plane 
motion  between  the  sphere  and  the  plane,  as  shown  in  Figure 
(3-12),  where  the  direction  of  the  sliding  force  Fs  and  the 
direction  of  the  friction  force  Ff  are  collinear  but  op- 
posite in  sense.  If  the  deviation  angle  6 is  smaller  than 
the  friction  angle  cp,  then  the  mechanism  will  not  move  at 
all.  Thus,  the  condition  that  the  sliding  force  must  over- 
come the  friction  force  should  be  satisfied  if  a continuous 
motion  is  desired.  It  is  also  observed  from  Figure  (3-12) 
that  a very  small  deviation  angle  may  not  be  desired  al- 
though the  smaller  is  the  deviation  angle,  the  better  is  the 
transmissibility  of  the  mechanisms. 

To  establish  the  general  condition,  considering  Figure 
(3-13),  the  sphere  and  the  plane  are  now  assumed  moving  in 
the  space.  The  x coordinate  is  assumed  to  be  in  the  direc- 
tion of  the  friction  force  and  the  y coordinate  is  in  the 
direction  of  the  normal  force,  as  shown  in  the  figure.  The 
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Figure  3-13  Friction  effects  in  spatial  motion. 
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absolute  magnitude  of  the  component  of  the  sliding  force  in 
the  direction  of  the  relative  motion  may  be  expressed  in 
terms  of  the  input  effort  force,  as 

I — sx I = | Ff | | sinScosaj  (3-42 

where  a is  the  angle  between  the  sliding  force  Fc„  and  the 
direction  of  the  relative  sliding  velocity  Vs.  The  mag- 
nitude of  the  friction  force  may  also  be  written 

Ilf  I = u|Fi | | cos6 | (3-43) 

where  u is  the  static  friction  coefficient.  Since  we  must 
have  the  condition  of  |FSX|  > |Ff|  to  satisfy  for  the  re- 
quirement for  continuous  motion  of  the  mechanism,  the  fol- 
lowing expression  is  given 

|cosa|  > u|cot6|  (3-44) 

In  the  other  words,  the  above  equation  may  also  be  employed 
to  restrict  the  allowable  friction  coefficient  between  the 
two  elements  in  contact,  i.e., 

U < |cosatan6|  (3-45) 

3,4  Mechanical  Advantage  and  Instantaneous  Efficiency 

The  friction  forces  will  cause  loss  of  energy,  which 
can  be  observed  directly  from  the  virtual  work  equation. 

For  example,  considering  the  mechanisms  containing  a one- 
point  contact  higher  pair,  the  virtual  work  equation  for 
these  mechanisms  can  be  expressed 

Md  • + •_Vd  + Ff  • Vs  + Mr  • «r  + Fr  • Vr  = 0 

(3-46) 

where  the  subscript  d and  r denotes  the  driving  link  and  the 
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resistant  link,  respectively;  m and  F denotes  the  torque  and 
force,  respectively;  Ff  is  the  sliding  friction  force;  V0  is 
the  relative  sliding  velocity  at  the  contact  point.  The 
first  two  terms  of  Equation  (3-46)  are  the  input  power 
(I.P.),  the  last  two  terms  is  the  output  power  (O.P.),  and 
the  power  loss  (P.L.)  is  then  given  by  the  third  term. 

Since  the  direction  of  the  friction  force  is  opposite  to  the 
direction  of  the  relative  sliding  velocity,  the  power  loss 
can  be  simplified  as 


P.L.  ~ | Ff | | Vs | (3-47) 

Observing  from  equation  (3-47),  the  power  loss  is  propor- 
tional to  the  magnitude  of  the  normal  force  and  the  mag- 
nitude of  the  relative  sliding  velocity.  The  power  ratio 

may  be  defined  as  the  ratio  of  the  output  power  to  the  input 
power,  as 


P.R. 


O.P. 


I.P. 


1 


P.L. 


I.P. 


(3-48) 


The  power  ratio  is  actually  the  instantaneous  efficiency 
(Alizade  and  Sandor,  1985).  The  instantaneous  efficiency 
may  be  defined  as 


where  1^*1  is  the  absolute  value  of  the  input  moment  re- 
quired without  friction  in  the  higher  pair  to  overcome  the 
output  resistance,  and  | Mjj | is  the  absolute  value  of  the  in- 
put moment  required  in  the  higher  pair  when  there  is  fric- 


tion. 
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The  mechanical  advantage  (M.A.)  is  one  of  the  criteria 
for  the  ability  of  a mechanism  to  transmit  torque  or  force 
and  is  defined  as  the  ratio  of  the  magnitudes  of  force  out 
over  force  in  (Erdman  & Sandor,  1984,  pp.  110  - 113): 


M.A.  = 


!out 


(3-50) 


m 


The  mechanical  advantage  may  also  be  investigated  on  the  in- 
formation of  input  and  output  displacements  (Midha,  etl., 
1984).  To  investigate  the  mechanical  advantage,  we  may  as- 
sume an  input  and  an  output  force  applied  at  a unit  normal 
distance  from  the  input  and  output  link  pivots,  respec- 
tively. if  the  output  joint  is  a cylindric  joint,  the  force 
output  has  two  components,  one  for  the  translational  and  one 
for  the  rotational  motion,  respectively.  Then,  we  may 
define  two  mechanical  advantages  for  each  of  the  force  out- 
put components,  i.e.. 


M. A. t = 


and 


M.A, 


FOUtfc 


m 


F 

rout 


m 


(3-51) 


(3-52) 


where  Foutt  and  Foutr  are  the  two  force  components  for  the 
translational  and  rotational  motion,  respectively.  These 
two  mechanical  advantages  must  satisfy  the  relationship  of 
M.A.  = (M.A.t2  + M. A . r 2 ) 1/2 


(3-53) 
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3 . 5 Acceleration  Ratio 

The  output  acceleration  ratio  can  be  defined  as  the 
ratio  of  the  maximum  and  minimum  accelerations  of  the  output 
link.  If  the  output  link  is  in  oscillation,  then  the  maxi- 
mum and  minimum  accelerations  happen  when  the  output  link  is 
at  the  limit  positions.  Thus,  we  have 

a = l^max  / ^minl  (3-54) 

where  cpmax  and  cpmin  are  the  maximum  and  minimum  output  an- 
gular or  linear  acceleration.  By  minimizing  the  quantity  jl 
a| , the  acceleration  can  be  forced  to  be  symmetrically 
distributed  in  the  positive  and  negative  directions.  To 
minimize  the  maximum  magnitude  of  the  acceleration,  another 
criterion  is  needed.  If  the  input  joint  is  a revolute  joint 
and  the  input  link  has  complete  rotation  and  constant  an- 
gular velocity,  such  other  criterion  can  be  defined  as  the 
ratio  between  the  maximum  magnitude  of  the  output  accelera- 
tion and  the  square  of  the  input  angular  velocity.  Thus, 

a*  = 1^1 max  / (3-55) 

When  the  constant  input  angular  velocity  is  unity,  the  ratio 
is  simply  the  maximum  output  acceleration  itself.  The  other 
way  to  minimize  the  maximum  magnitude  of  the  maximum  ac- 
celeration is  to  minimize  the  deviation  between  the  maximum 
and  minimum  accelerations.  Thus,  we  may  define  the  ac- 
celeration deviation,  as 

a = l(cpmax  " ^min ) /(Pmax  I (3-56) 

Equations  (3-54)  and  (3-56)  may  be  also  applied  to  output 
linear  acceleration. 


CHAPTER  4 

STATIC  FORCE  AND  TORQUE  ANALYSIS 

The  static  force  and  torque  analysis  is  carried  out  by 
vector  methods  for  three  numerical  examples.  The  first  one 
is  an  R-Sp-R  mechanism,  which  has  one  contact  force  in  the 
model  of  the  Sp  pair  and  is  statically  determinate.  The 
second  one  is  an  R-Cp-C  mechanism,  which  has  two  contact 
forces  in  the  model  of  the  Cp  pair.  The  R-Cp-C  mechanism  is 
a statically  indeterminate  mechanism,  therefore,  the  output 
force  and  moment  are  specified  instead  of  the  input  moment. 
The  last  one  is  an  R-PL-C-R  mechanism,  which  has  three  con- 
tact forces  in  the  model  of  the  PL  pair  and  is  also  stati- 
cally determinate.  The  static  force  and  torque  analysis  for 
other  mechanisms  can  be  followed  in  a similar  manner. 

4 • 1 R-Sp-R  Mechanisms ( Figure  (4-1)) 

In  the  case  of  statically  determinate  three-link 
mechanisms  such  as  the  R-Sp-R  mechanism  with  one  point  con- 
tact between  the  higher  pair  elements,  which  is  shown  in 
Figure  (4-1),  the  moment  equilibrium  equation  for  the  input 
link  can  be  written  as 

(Md  + P2  x R32 ) • u2  = 0 (4-i) 

where  is  the  driving  torque;  R32  is  the  reaction  force  at 
the  contact  point;  P2  is  the  vector  from  the  input  R joint 
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Figure 


1 Free  body  diagram  of  R-Sp-R  mechanism. 
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to  the  contact  point;  and  u2  is  the  unit  vector  of  the  input 
R joint  axis.  The  reaction  force  consists  of  the  trans- 
mitted normal  force,  Fn32,  and  the  sliding  friction  force. 

If  32*  This  can  be  expressed,  as 

—32  = In32  + If 3 2 

= Hn32l (sign* 1(fn32)fn32  + Uff32)  (4-2) 

where 


— n32 

= - a 

If  32 

= - Is 

^s 

= 

Ys/|YS| 

w 

>1 

= 

Yp'p  (YP 

'P 

= V - Yp 

V 

! t 

= £3  x I3 

a 

>1 

= 

«2  x 12 

a)a 


The  magnitude  of  the  normal  force  is  given  by  substituting 
Equation  (4-2)  into  Equation  (4-1),  i.e., 

- IVU 


Fn  = 


or 


12  x ( sign( fn32 ) fn32  + Uff32)  * u2 
- Md 


(4-3) 


(4-4) 


12  x r32  ‘ H2 

where  r32  = sign(  fn32  )coscpfn32  + sincpff32  and  R is  the  mag- 


1.  The  direction  of  the  normal  force  applied  from  the 
cylinder  to  the  plane  is  initially  defined  the  same  as  the 
normal  vector  of  the  contact  plane,  fn23  = a,  i.e,  fn35  = 
a.  Then,  it  is  verified  by  the  sign  of  the  absolute  mag- 
nitude of  the  normal  force,  | Fn | . if  the  resulting  I F_ I is 
positive  then  it  agrees  with  the  assumption.  A negative 
I In  I means  that  fn32  = a. 
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nitude  of  the  reaction  force,  R32.  The  above  equation  is 
then  employed  to  check  the  sign  of  fn32.  The  moment  equi- 
librium equation  for  the  output  link  can  be  written 

(Mr  + P3  x R23)  • u3  = 0 ( 4 — 5 ) 

where  Mr  is  the  output  resistant  torque,  P3  is  the  distance 
vector  from  the  output  joint  to  the  contact  point.  Thus, 
the  resistant  torque  is  given  by  substituting  Equation  (4-4) 
into  the  above  equation,  as 


- Md(P3  x r32  • u3) 


To  derive  the  instantaneous  efficiency  to  show  the  fric- 
tional effects  explicitly,  we  may  write  the  virtual  work 
equation  for  the  mechanism,  as 

^d  * 9.2  + £f  • Ys  + Mr  * ^3  = 0 (, 

— f is  the  friction  force  where  Ff  = -u|Fn|vs.  The  power 
loss  is 


= -sin<p|R|  |VS| 

Now,  the  instantaneous  efficiency  can  be  written  in  terms  of 
the  power  loss  and  the  input  power,  as 


I.E.  = 1 


where 


P.L. 


I.P. 


(4-9) 


I • P • - • fi2 

Substituting  Equation  (4-8)  into  the  above  equation  gives 
the  instantaneous  efficiency  in  terms  of  the  geometrical 


parameters,  as 
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I.E.  = 


slncp  | Vs  | 


(4-10) 


—2  x r32  • Q2 

To  investigate  the  mechanical  advantage,  we  may  assume 
a force  input  and  a force  output  applied  at  the  input  and 
output  link,  respectively.  The  directions  of  the  force  in- 
put and  force  output  are  assumed  to  be  the  same  as  the 
directions  of  the  absolute  velocities  of  the  contact  point  P 
as  a point  of  the  input  link  and  the  output  link,  respec- 
tively. Then,  the  force  input  and  force  output  may  be  writ- 
ten in  terms  of  the  magnitudes  of  the  input  torque  and 
resistant  torque,  respectively,  as 

iMdl 


I — i I = 


Hi  I 


(4-11) 


and 


'-OI  I Pol  (4-12> 

where  |FjJ  and  |F0|  are  the  absolute  magnitude  of  the  force 
input  and  the  force  output,  respectively;  |PjJ  and  |P  | are 
the  absolute  magnitude  of  shortest  distance  from  the  input 
and  output  joint  axes  to  the  force  input  and  output,  respec- 
tively. Thus,  the  mechanical  advantage  for  this  particular 
mechanism  can  be  written  in  terms  of  the  input  and  output 
moment , as 


M.  A. 


I Mr | | P j 1 
l*Sdl  Hoi 


(4-13) 


or  can  be  written  in  terms  of  the  instantaneous  efficiency, 
as 
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M.A.  = 


I ^2 II  Pi  I 


I .E. 


(4-14) 


I —3  I I Eo  I 

If  more  mechanical  advantage  is  desired,  one  may  simply  in- 
crease |Pj_|,  i.e.,  the  distance  from  the  input  joint  axis  to 
the  input  force;  or  decrease  |PQ|,  i.e,  the  distance  from 
the  output  joint  axis  to  the  output  force. 


4 • 2 R-Cp-C  Mechanisms ( Figure  (4-2)) 

For  the  R-Cp-C  mechanisms,  which  is  a statically  in- 
determinate mechanism,  we  may  model  it  as  a function  genera- 
tor, i.e.,  the  output  link  is  oscillating  in  a limited 
range,  and  specify  the  output  resistant  force  and  torque. 

In  order  to  do  so,  the  output  cylinder  joint  is  coupled  with 
a translation  compression  or  tension)  spring  and  a torsion 
spring,  as  shown  in  Figure  (4-2).  Assuming  the  cylinder 
joint  is  free  of  friction,  then  the  resistant  force  and  mo- 
ment are  simply  proportional  to  the  linear  and  angular  dis- 
placements, respectively,  as 

-r  = kt(_S3)  (4-15) 

and 


where  kt  and  kr  are  the  translation  and  torsion  spring  con- 
stants, respectively;  S3  and  03  are  the  linear  and  angular 
displacements  of  the  output  joint,  respectively. 

Now,  if  the  output  motion  is  known,  the  resistant  force 
and  moment  would  be  given,  we  may  write  the  force  equi- 
librium for  the  output  link  as  follows: 
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Figure  4-2  Free-body  diagram  of  R-Cp-R  mechanism. 
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Md  t(clt>2  “ C2b^)Fr  + + C2a]_)Mr]/D  (4-22) 

where 

Ci  = -2i  x (sign(fni23)coscPfni23  + sincpffi23)  • u2; 
i = 1,  2 

Equation  (4-22)  shows  that  the  driving  torque  is  determined 
by  two  parts,  one  is  by  the  resistant  force  and  the  other  is 
by  the  resistant  moment.  To  illustrate  the  effects  of  the 

Coulomb  friction,  we  may  write  the  virtual  work  equation  for 
the  mechanism,  as 

Md  • - sincp(  iRxl  | Vsl  | + | R2  | | Ys2  I ) + Mr  • + 

-r  ' -3  = 0 (4-23) 

where  |Vsl|  and  |Vs2|  are  the  absolute  magnitude  of  the 
relative  sliding  velocities  at  the  two  contact  points  of  the 
Cp  pair,  respectively,  and  and  S3  are  the  angular 
velocities  of  the  input  and  the  output  links.  Thus,  the 
power  loss  is  explicitly  expressed,  as 

P.L.  = sincp(  | | | Vsl  | + |R2  I | Vs2  | ) (4-24) 

Substituting  Equations  (4-19)  and  (4-20)  into  the  above 
equation  gives  the  P.L.  in  terms  of  the  magnitude  of  the 
resistant  force  and  moment,  as 

P.L.  = sincp[(-b2|Vsl|  + b2  |Vs2 1 )Fr  + [a2\Vsl\  - 

al IXs2 I )^r J /D  (4-25) 

The  output  power  is  consisted  of  two  components  and  the  in- 
stantaneous efficiency  may  be  expressed  as  the  summation  of 
two  instantaneous  efficiencies  for  output  force  and  output 
moment,  respectively,  as 

I.E.  = I.E.p  + I . E . jyj 


(4-26) 
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where 


I . E . p 


and 


0 . P . p 

1. P. 


Frv3 

Mdfi2 


I-P-  MdS22 

Here,  instead  of  using  the  absolute  magnitudes,  the  signed 
scalar  are  used,  since  the  output  force  and  moment  may  con' 
tribute  either  positive  or  negative  work.  To  express  the 
total  instantaneous  efficiency  in  terms  of  the  total  power 
loss,  one  may  write 


I.E. 


where 


O.P. 

Mrn3  + Frv3 

P.L. 

I.P. 

Mdfi2 

I.P. 

F*L*  (AlFr  + A2Mr)sincp 

I-P.  B1Fr  + B2Mr 

A1  = ~b2 I Ysi | + b2|Vs2| 

A2  = a2 I Ysl I - allYs2| 

B1  = ^2 1 — pi  * r132|  - b1|vp2  • r232| 

b2  = -a2 I — pi  * £132  I + allYp2  • r232 | 

The  mechanical  advantage  may  also  be  expressed  in  two 
terms  but  not  by  the  way  of  superposition  as  the  power 
ratio.  This  is  because  the  force  output  consists  of  two 
components:  one  is  favoring  in  the  translational  motion  and 
the  other  one  is  favoring  the  rotational  motion.  The  mag- 
nitude of  the  resultant  force  is  the  square  root  of  the  sum 
of  the  squares  of  the  magnitudes  of  force 


components . Thus , 
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the  total  mechanical  advantage  may  be  written  in  terms  of 
the  mechanical  advantages  for  force  output  produced  by  out 
put  force  and  output  moment,  respectively,  as 

M.A.  = (M. A. t2  + M.  A. r2 ) !/2  (4_2l 

where 


M.A.t 


1 — r I I £i  I 
iMdl 


m.  a.  r = = JMlzil 

I— i I |Mdl|Pol 

The  absolute  magnitudes  | Fot | and  | For | are  the  force  output 
components  for  translational  motion  and  rotational  motion, 
respectively;  iPjJ  and  |P0|  are  the  distances  from  the  joint 
axes  to  the  locations  of  the  force  input  and  force  output. 
Unlike  the  previous  case,  the  power  ratio  and  mechanical  ad- 
vantages for  the  R-Cp-C  mechanism  can  not  be  expressed  ex- 
plicitly in  terms  of  velocities. 


4 • 3 RzPt.~C~R  Mechanism)  Figure  4-3) 

Considering  an  r-pl-c-r  mechanism,  as  shown  in  Figure 
(4-3),  the  moment  equilibrium  for  the  input  link  can  be  ex- 
pressed 

Md  + (£21  x R132  + P22  x R232  + p23  x r332)  . u2  = 0 

(4-29) 

where  P2i,  i = 1,  2,  3,  are  the  vectors  from  u2  to  the  three 
contact  points  (two  in  the  outter  plane  element  and  the  one 
in  the  inner  plane  element),  respectively.  Recognize  that 
if  there  is  no  force 


or  moment  transmitted  in  the  direction 
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*32b4> 

-32b  f 


Figure  4-1  Free-body  diagram  of  R-P 


L~C-R  mechanism. 
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°f  -3  (i*e*'  there  is  no  friction),  then  from  the  force  and 
moment  equilibrium  equations  for  the  coupler  link (link  3)  we 
have 

^ — 123  + 1223  + 1323  ) * £3  = 0 (4-30) 

(— 31  x —123  + 132  x 1223  + 133  x 1323 ) * I3  = 0 (4-31) 

Then,  the  maqnitude  of  the  reaction  forces  can  be  obtained 
by  solvinq  Equations  (4-29)  to  (4-31)  simultaneously.  Usinq 
the  above  three  equations  we  obtain  three  equations  linear 
in  the  three  absolute  maqnitudes  of  the  normal  forces,  as 

alllll  + ll2 I + c 1 1 13 | = -Md  (4-32) 

aillll  + ll2 I + ci II3 | =0  i = 2,  3 (4-33) 

where 

al  = 121  x Ii32  * £2 

bl  = 122  x 1232  • 112 
C1  = 123  x £332  * £2 
a2  = £123  ' £3 
b2  = £223  * £3 
c2  = £223  * £3 
a3  = 131  x £123  ‘ £3 
b3  = 132  x £223  ‘ £3 
c3  = 133  x £323  * £3 

— i32  ~ sign(fni32  )cos<Pf.ni32  + sincpffj_32  = - £123; 
i = 1,  2,  3 

Applying  Cramer's  rule,  the  reaction  forces  can  be  expressed 
llll  = MdDl/D  (4_34) 

1-2  I = MdD2/D  . 


and 


|R3|  = MdD3/D 
respectively;  where 
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(4-36) 


D 


al  bl 
a2  b2 
a3  b3 


C1 

c2 

c3 


D1  = b2c3  " b3c2 
D2  = a3c2  - a2c3 
d3  = a2b3  “ a3b2 


Equations  (4-34)  to  (4-36)  are  employed  to  check  to  the  sign 
of  the  normal  forces.  The  resistant  moment  can  be  obtained 
by  solving  the  moment  equilibrium  equation  for  the  dyad  con- 
sisting of  the  coupler  link  and  the  output  link,  as 

^—41  x —123  + £42  x ^223  + £43  x ^323)  ' u4  + 

-r  * -4  = 0 (4-37) 

Thus,  by  substituting  Equations  (4-34)  to  (4-36)  into  the 
above  equation,  the  resistant  moment  can  be  expressed  in 
terms  of  the  input  moment,  as 

Mr  = (£41  x £.123  )D1  + (P42  x £223 )D2  + 

(—43  x £323  ^ D3  ^Md  * vx4  > /D  (4-38) 


To  derive  the  instantaneous  efficiency,  the  virtual 
work  equation  for  the  mechanism  may  be  written 

Md  • - sincp(  | Rx  | | Vsl  | + |R2  I | Vs2  | + | R3  | |Vs3  | ) + 

— r * ^4  = 0 ( 4_ 

Thus,  the  instantaneous  efficiency  for  this  particular 
mechanism  may  be  expressed  in  terms  of  the  power  loss,  as 


O.P. 

I«rl 

l«4  | 

P.L. 

I.p. 

l«dl 

|n2l 

“ — 1 - 

I.P. 

where  the  ratio  of  power  loss  to  the  input  power  is 


(4-40) 
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P.L.  ^ UDlllYsll  + I ^-*2  I I — s2  I + Id3  I I Ys3  1 ) sin<P 
I.P.  |dq2| 

(4-41) 

The  mechanical  advantage  may  also  be  written  in  terms  of  the 
instantaneous  efficiency,  as 

Hoi  |n2ll£il 

M.A.  = = I.E.  (4-42) 

I Hi I I «4 I I Pol 

where  and  the  input  force  and  output  force,  respec- 
tively. 


CHAPTER  5 

DESIGN  METHOD  AND  FORMULATION 

5 . 1 Optimization  Strategy  in  Mechanisms  Design 
The  application  of  optimization  in  mechanism  design  may 
be  classified  into  two  groups,  one  is  in  precision  point 
synthesis  and  the  other  one  is  the  mathematical  programming 
problem  (Fox  & Gupta,  1973).  This  work  naturally  falls  into 
the  second  group.  The  optimization  methods  may  also  be 
classified  into  three  groups:  objective  function  without 
constraints,  objective  function  with  constraints  using  un- 
constrained optimization  methods,  and  the  direct  solution  of 
constrained  problems.  The  second  group  may  be  called  in- 
direct methods  for  constrained  problems  or  sequential  un- 
constrained minimization  techniques  (SUMT)  (Fiacco  & McCor- 
mick, 1968).  The  indirect  methods  handle  the  constraints 
implicitly  so  that  the  detailed  structure  of  the  constraint 
equations  does  not  need  to  be  considered.  The  third  group 
may  be  called  direct  solution  methods  of  constrained 
problems.  The  direct  methods  take  the  constraint  equations 
directly  into  account.  Most  application  of  optimal  en- 
gineering design  involve  nonlinear  equation  systems  with 
both  linear  and  nonlinear  constraint  equations.  Thus, 
generally  the  second  or  the  third  method  is  needed.  The 
comparison  of  these  methods  for  constrained  problems  has 
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been  studied  by  several  publications,  such  as  the  one  by 
Eason  and  Fenton  (1974)  and  more  recently  by  Sandgren  and 
Ragsdell  (1980).  A summary  for  these  and  other  comparison 
is  seen  in  Reklaitis,  Ravindran,  and  Ragsdell,  1983.  Hock 
and  Schittkowski  (1983)  also  compared  27  nonlinear  program- 
ming codes.  For  those  problems  where  the  gradients  are 
available  for  the  constraint  equations,  the  sequential  quad- 
ratic programming  method  and  the  generalized  reduced 
gradient  method  may  give  the  best  approach (Hock  and  Schit- 
tkowski, 1983).  Otherwise,  the  indirect  methods  will  be  a 
better  choice.  Although,  the  gradients  can  be  obtained  by 
numerical  methods,  it  does  not  necessary  give  a better 
result,  as  indicated  by  Eason  and  Fenton  (1974). 

The  general  optimal  design  problem  may  have  the  form 
(Rao,  1978) 


Minimize  U(X) 

Subject  to 


— ^X1 ' x2 ' • • • f xn^ 


T 


Gj(X)  >0  j = [1,  J] 

Hk(X)  =0  k = [J+l,  K] 

where  U(X)  is  the  objective  function,  X is  the  design  vec- 
tor, Gj(X)  are  the  inequality  constraints,  and  Hk(X)  are  the 
equality  constraints.  To  use  the  penalty  function  with  an 
unconstrained  optimization  method,  the  original  objective 
function  and  the  constraints  are  converted  into  the  form 

tjm(X,  R)  = U(X)  + M(R,  G ( X ) , H(X)  ) (5-D 

where  R is  a set  of  penalty  parameters  and  M is  a function 
of  R and  the  constraint  equations.  Now,  the  constrained  op- 
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tiraization  problem  is  treated  as  an  unconstrained  optimiza- 
tion problem  with  the  objective  function  uM. 

Powell's  method  using  conjugate  directions  (1964)  for 
the  unconstrained  optimization  with  the  method  of  penalty 
multipliers  (Powell,  1969;  Fletcher,  1974)  for  the  penalty 
function,  provides  an  efficient  optimization  scheme  for  con- 
strained optimization  without  the  need  of  derivatives.  The 
codes  form  OPTIVAR  (Siddall,  1982)  are  adopted  here  due  to 
its  efficiency(Reklaitis,  Ravindran,  and  Ragsdell,  1983)  and 
easy  implementation  with  different  penalty  functions.  The 
objective  function  illustrated  by  Siddall  has  the  form 

1 1 

UM  - U(X)  + - E r j <Gj  - gj>2  + _ 2 rk(Hk  _ hk ) 2 (5_2) 

^3  2 k 

where  rj  and  rk  are  the  multipliers  and  may  vary  for  succes- 
sive optimization  runs,  which  are  recommended  to  be  in- 
creased by  an  arbitrary  factor  of  10.  The  <*>  notation  has 
the  meaning 

* if  * < 0 

< *>  = 

0 if  * > o (5"3) 

The  other  two  parameters  gj  and  hk  are  updated  by  the  fol- 
lowing equations  between  successive  iteration  i and  i+1,  as 

9j(l+1)  = gj(i)  - min  (G-j  , gj<i>)  (5-4) 

hk(1+1)  = hk(i>  - HkU)  (5.5) 

The  first  run  uses  gj  and  hk  equal  to  zero.  Figure  (5-1) 
illustrates  a flow  chart  adapted  from  the  algorithm  listed 
by  Siddall  for  Powell's  strategy. 
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Quantities  without  supscripts  are  current  values. 


Figure  5 1 Flow  chart  of  Powell's  algorithm  for  method 
multipliers  (adapted  from  Siddall,  1982). 


of 
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5 • 2 Problem  Formulation 

Involving  the  transmission  criterion  in  optimal  design 
of  mechanisms  has  been  seen  in  numerous  publications.  Some 
of  them  work  on  finding  the  optimal  transmission  angle  or 
deviation  angles  for  a certain  mechanism,  some  others  take 
the  transmission  angle  as  one  of  the  design  criteria.  Here, 
the  goal  is  to  find  the  best  transmission  characteristics 
which  includes  not  only  minimizing  the  maximum  deviation 
angles  but  also  minimizing  the  peak  acceleration  ratio. 
Having  the  transmission  criteria  defined  in  the  previous 
chapters,  we  need  to  formulate  the  design  problem  for  each 
particular  mechanism.  Before  doing  so,  some  basic  con- 
straints need  to  be  clarified  first. 

5.2.1  Rotatability  of  mechanisms 

To  investigate  the  transmission  characteristics,  a 
fully  rotatable  input  link  is  desired,  therefore  the 
rotatability  of  the  mechanism  must  be  considered.  The  term 
"rotatability"  applies  to  not  only  the  input  link  but  also 
the  output  link.  Thus,  multiple  conditions  of  rotatability 
should  be  established  to  cover  different  Grashof-type 
linkages  namely  crank-rocker  or  crank-crank. 

Generally,  the  rotatability  for  the  mechanisms  having 
quadratic  input-output  displacement  equation  may  be  deter- 
mined by  the  natures  of  the  roots  of  equation's  dis- 
criminant, as  seen  in  Nolle 's  work  (1969)  for  RSSR 
mechanism.  Alizade  and  Sandor  (1985)  also 


apply  the  ine- 
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quality  constraints  derived  from  the  coefficients  of  the 
determinant  to  the  design  of  RSSR  and  RSSP  mechanisms  with 
complete  crank  rotation.  Reinholtz  (1983)  uses  a similar 
approach  to  determine  the  input  crank  rotatability  for  the 
RSSR  mechanism.  Recently,  Angeles  and  Bernier  (1986)  use 
the  discriminant  of  the  input-output  equations  by  geometric 
approach  to  the  analyses  of  the  mobility  of  planar  and  spa- 
tial four-bar  mechanisms.  The  determining  of  rotatability 
for  the  R-Sp-R  mechanism  (see  Figure  5-2)  will  need  to  util- 
ize a similar  approach. 

5. 2. 1.1  Rotatability  of  R-Sp-R  mechanism 

The  discriminant  a quadratic  displacement  equation  can 
be  rearranged  into  a quartic  equation  in  terms  of  either  the 
output  or  input  rotation  angle.  After  the  discriminant  of 
the  quartic  equation  is  obtained,  Sturm's  theory  may  be  ap- 
plied to  judge  the  roots  of  the  equation  (Cajori,  1943).  If 
the  discriminant  of  the  quartic  equation  has  four  real 
roots,  then  the  original  displacement  equation  must  have 
complex  roots  only,  i.e.,  the  input  link  must  have  no  limit 
positions,  which  also  means  the  input  link  has  a complete 
rotation.  On  the  other  hand,  if  the  discriminant  of  the 
quartic  equation  has  no  real  roots  at  all,  the  original  dis- 
placement equation  must  have  real  roots,  i.e.,  there  are 
limit  position,  for  the  input  link.  In  another  words,  the 
input  link  does  not  have  full  rotatability.  Similar  situa- 
tions are  applied  to  the  output-input  equation  to  determine 
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Figure  5-2  An  R-Sp-R  mechanism 
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the  rotatability  of  the  output  link. 

To  state  the  criteria  for  the  quartic  equation  having 
either  all  real  roots  or  all  imaginary  roots,  some  brief  ex- 
pressions from  Cajori's  work  (1943,  Sec.  50)  are  adopted  as 
follows : 

A general  quartic  equation  may  have  the  form 
b0x4  + 4b1x3  + 6b2x2  + 4b3x  + b4  = 0 (5-6) 

which  can  be  transformed  into  a new  equation,  deprived  of 
its  second  term  as  follows: 

y4  + 6Hy2  + 4 Gy  + b02I  - 3H2  = 0 (5-7) 

where 

y = b0x  + bx 

H = b0b2  “ bl2 
G = bgb3  - 3bQb^b2  + 2b^3 

I = bob4  ” 4b^b3  + 3b22 

The  nature  of  the  roots  of  the  general  quartic  equation  can 
then  be  stated  as  follows: 

a.  all  roots  real,  if 

(I3  - 27 J2 ) > 0, 

(3b0J  - 2HI)  > 0, 
and 
H < 0; 

b.  all  roots  imaginary,  if 

(I3  - 27 J2 ) > 0, 
and  if 
H > 0 


or  else 
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(3b0J  - 2HI ) < 0; 

where 

J = (b02HI  - G2  - 4H3)/b03. 

To  apply  the  above  two  conditions  in  mechanism  design, 
we  may  affirm  the  following  conditions: 

1.  if  the  determinant  of  the  quadratic  output-input 

equation  of  a mechanism  satisfies  the  first  condi- 
tion, then  the  output  link  is  a rocker,  i.e.,  it 
is  oscillating  in  a limited  range; 

2.  if  the  determinant  of  the  quadratic  output-input 

equation  of  a mechanism  satisfies  the  second  con- 
dition, then  the  output  link  is  a crank,  i.e.,  it 
can  rotate  360  degrees; 

3 . if  the  determinant  of  the  quadratic  input-output 

equation  of  a mechanism  satisfies  the  second  con- 
dition, then  the  input  link  is  a crank,  i.e.,  it 
has  full  rotatability . 

Generally,  these  conditions  can  be  applied  to  the  four-link 
spatial  mechanisms  with  lower  pairs  only  and  the  three-link 
spatial  mechanisms  with  one  higher  pair. 

5. 2. 1.2  Rotatability  of  R-Cp-C  and  R-PL-C-R  mechanisms 

It  is  noted  that  the  R-Cp-C  and  R-PL-C-R  mechanisms,  as 
far  as  rotations  are  concerned,  are  equivalent  to  a spheri- 
cal RRRR  four-link  mechanism  as  shown  in  Figure  5-3.  The 
parameters  used  to  specify  it  are  as  follows: 

£ the  unit  vector  along  the  ith  pair  axis 
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Figure  53  A spherical  RRRR  four-link  mechanism. 


Wi 
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a.j_j  the  unit  vector  along  the  common  perpendicular  be- 
tween Sj_  and  s j 

a^j  the  positive  angle  from  sj_  to  sj  measured  clock- 
wise looking  in  the  direction  on  j , i.e.,  in  the 
right-hand  direction  about  a^j 
©i  the  positive  angle  from  a-jj  to  a.j^  measured  clock- 
wise looking  in  the  direction  of  Sj 
In  the  equivalent  mechanism,  the  a.jj's  become  arcs  of  great 
circles  physically  representing  the  link  lengths.  Since  an- 
gular relationships  are  maintained,  the  spherical  mechanism 
has  the  same  Grashof  type  as  the  spatial  mechanism. 

The  input-output  equation  for  the  4R  spherical 
mechanism  is  a cosine  law  for  a spherical  quadrilateral 
(Duffy,  1980).  If  all  the  link  lengths  lie  within  certain 
ranges,  the  Grashof  rule  for  the  planar  four-link  four- 
revolute  mechanism  can  be  directly  applied  to  the  spherical 
four-link  four-revolute  mechanism  using  the  arc  angles  be- 
tween joints  as  link  lengths.  Gilmartin  and  Duffy  (1972) 
derived  the  criteria  for  determining  types  of  the  spherical 
four-link  four-revolute  mechanism  with  general  proportions. 
The  criteria  for  the  crank-rocker  type  mechanism  are  stated 
as  three  inequalities: 


a23 

+ 

a34 

> 

a12 

+ 

a41 

(5-8) 

a23 

+ 

a41 

> 

a12 

+ 

a34 

(5-9) 

a3  4 

+ 

a41 

> 

a12 

+ 

a23 

(5-10) 

The  criteria  for  the  crank-crank  type  mechanism  are  stated 
as  another  three  inequalities: 
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a12 

+ 

a23 

> 

a34 

+ 

a41 

(5-11) 

a23 

+ 

a34 

> 

a23 

+ 

a41 

(5-12) 

a23 

+ 

a34 

> 

a12 

+ 

a34 

(5-13) 

To  utilize  these  inequalities,  we  must  first  observe  that 
the  R-Cp-C  and  R-PL-C-R  mechanisms  are  equivalent  to  the 
spherical  four-link  four-revolute  mechanism.  This  can  be 
observed  by  modeling  all  the  pairs  with  lower  pairs  only. 
Figure  5-4  shows  an  R-Cp-C  mechanism.  The  Cp  pair  can  be 
modeled  by  the  C-P-R  dyad.  Figure  5-5  shows  an  R-CPR-C 
mechanism  which  is  equivalent  to  the  R-Cp-C  mechanism.  In 
Figure  5-5,  the  axes  of  the  turning  pairs  are  s^,  S2,  S4, 
and  s5.  These  axes  are  corresponding  to  u2,  £*/  a,  and  u.3, 
respectively,  as  shown  in  Figure  5-4.  The  input-output 
equation  for  the  angular  motions  is  related  to  the  input  and 
output  angles  and  the  angles  between  the  axes  of  the 
revolute  pairs  only.  This  can  also  be  observed  from  the 
displacement  equation  (E-3)  in  Appendix  E.  To  obtain  a 
crank-rocker  type  mechanism,  equations  (5-8)  to  (5-10)  are 
applied.  With  proper  substitutions  of  the  corresponding 
angles  between  the  axes , the  inequality  constraints  can  be 
expressed  as: 

a24  + a45  " ( a12  + a51 ) > 0 (5-14) 

a24  + a51  " ( a12  + a45 ) > 0 (5-15) 

a45  + a51  " ( a12  + a24)  > 0 (5-16) 

Notice  that  in  Figure  5-5  the  angle  between  s2  and  £4,  a24, 

is  the  same  angle  between  c and  a in  Figure  5-4,  which  is 
n/2.  The  other  angles  can  be  expressed  as: 
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Figure  5-4  An  R-Cp-C  mechanism. 
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Figure  5-5  An  R-CPR-C  mechanism. 
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a-12  = cos-l 

£2 

* £ 

(5-17) 

a45  = 

a • 

£3 

(5-18) 

a51  = cos--*- 

£3 

* £2 

(5-19) 

Similarly,  an  R-PL-C-R  mechanism  (shown  in  Figure  5-6) 
can  be  modeled  as  an  R-PRP-C-R  mechanism  (shown  in  Figure 
5-7).  The  axes  of  the  turning  pairs  in  Figure  5-7  are  slf 
£3  • £5'  and  £6'  which  are  corresponding  to  u2,  a,  u3 , and 
u4,  respectively.  To  obtain  a crank-crank  type  mechanism, 
equations  (5-11)  to  (5-13)  are  applied.  With  proper  sub- 
stitutions of  the  corresponding  a angles,  the  inequality 
constraints  can  be  written  as: 

a13  + a35  ~ ( a56  + a6l)  > 0 (5-20) 

a13  + a56  " (a35  + a6l)  > 0 (5-21) 

a35  + a56  " ( a12  + a41 ) > 0 (5-22) 

where  the  a angles  can  be  obtained  from  observing  Figure 


as 

a13  = cos-1 

£2 

• a 

(5-23) 

0,35  = cos--'- 

a • 

£3 

(5-24) 

0.55  = cos-1 

£3 

* £4 

(5-25) 

a61  = 

£4 

* £2 

(5-26) 

5.2.2  Formulation  of  R-Sp-R  mechanism 

An  R-Sp-R  mechanism  is  shown  in  Figure  (5-2).  The  in- 
put joint  axis  lies  in  the  x coordinate  axis,  the  output 
joint  axis  is  a unit  distance  away  from  the  x coordinate 
axis  with  a skew  angle  about  the  y coordinate.  Thus  the 
common  normal  of  the  input  and  output  joint  axes  is 
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Figure  5-6 


An  R-Pl-c-R  mechanism. 
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a-34  s.4 


Figure  5-7  An  R-PRP-C-R  mechanism. 


Ill 


[0,1,0]-'-.  The  skew  angle  of  the  output  joint  axis  may  be 
assigned  as  a fixed  value.  The  position  vector  P locates 
the  initial  contact  point  and  the  shortest  vector  from  the 
input  joint  axis  to  the  contact  point  is  denoted  as  P2.  A 
coordinate  system  attached  to  the  contact  point  has  axis  a 
in  the  direction  of  the  normal  vector  of  the  plane  element 
of  the  Sp  pair,  axes  b and  c are  on  the  contact  surface. 

The  goal  is  to  design  a crank-crank  mechanism  with  favorable 
transmission  criteria  for  a specified  link  length  ratio  be- 
tween the  input  crank  and  the  ground  link.  The  state  vari- 
ables are 

M.2'  u 3 ' 22,  D3,  P2,  a,  b,  c 

where 

a • b = 0 ; 
c = a x b; 

and 

D2  = [0,  0,  0]T, 

D3  = [0,  1,  0]T 

are  the  vectors  from  the  origin  of  the  fixed  coordinate  sys- 
tem to  the  input  joint  and  output  joint,  respectively.  The 
design  variables  are  the  position  vector  P and  the  Euler 
angles  of  the  unit  vector  a.  The  use  of  Euler  angles  for 
the  in-plane  coordinate  axes  has  the  advantage  of  limiting 
the  ranges  of  the  magnitudes  of  the  design  variables.  Thus 
the  design  vector  may  be  expressed 

— - t ®x ' ®y'  ®z r Px'  Py'  PZ]T 
Then , we  have 


(5-27) 
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I.  f-Px'  Py  / Pz  ^ 
a i = [ Mq ] a 
ki  = [M0]b 

Ci  = ai  x bj_ 

where  [M0]  is  the  rotation  matrix  based  on  ©v,  ©^ ; ©v 

a.'  y r ca  r xf 

0y,  0Z  are  the  angles (at  the  i-th  iteration)  of  the  moving 
coordinate  system  measured  from  the  positive  branch  of  the 
fixed  coordinate  axes,  respectively;  and  a^,  bj_,  and  c^ 
designate  the  moving  coordinate  axes  at  the  i-th  iteration. 
Then,  the  rotation  matrix  for  the  moving  coordinate  system 
at  the  i-th  iteration  can  be  expressed 

1 

-C@yS0x  S0y 

-s0xs©ys©z+c©zc©x  -s0zc©y 
S0XS@yC0Z+S©ZC0X  CQZC0y 

(5-28) 

where  c©j  = cosQj  and  s©j  = sinQj ; 0j  = 0X,  0y,  0Z. 

5. 2. 2.1  Objective  function 

The  objective  function  may  be  based  on  the  deviation 
angles  and  the  maximum  and  minimum  acceleration  of  the  out- 
put joint,  as 

U(X)  = - k1|cos6F23 |min  - k2 lcos6FR|min  + k3aQ3  (5-29) 

where 

X = [0X,  0y,  0Z,  Px,  Py,  Pz]T. 

<*93  = I ®3  I max  / 022- 

©2  and  03  are  the  input  angular  velocity  and  the  output  an- 


CMq]  = 


C©yCQX 


C0xS©yS0z+C©zS0x 

-C0xS0yC0z+S©zS0x 
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gular  acceleration,  respectively.  The  coefficients  k-j_,  k2, 
and  k3  are  the  weighting  factors.  The  two  deviation  angles 
are  defined  in  Equations  (3-19)  and  (3-20). 

5. 2. 2. 2 Inequality  constraints 

The  R-Sp-R  mechanism  is  desired  to  be  a crank-crank 
type  mechanism,  so  that  the  second  and  third  conditions  of 
the  rotatability  are  to  be  applied.  Then,  for  the  input- 
output  and  output-input  equations  we  have  six  inequality 
constraint  equations,  i.e.. 


Gi(X)  = (Ii3  - 27Ji2)  > 0; 

(5-30) 

g2^)  = Hj_  > 0 or  else 

(5-31) 

G3(X)  = ObiQJi  - 2HiIi)  < 0; 

(5-32) 

G4(X)  = (IQ3  - 27Jq2)  > 0; 

(5-34) 

Gs(X)  = H0  > 0 or  else 

(5-35) 

G6(X)  = (3bo0Jo  - 2H0I0)  < 0; 

(5-36) 

where  the  subscripts  i and  o denote  that  the  coefficients 
are  derived  from  the  coefficients  of  the  input-output  and 
output- input  equations,  respectively.  The  coefficients  of 
the  quartic  equations  are  derived  in  Appendix  G. 

As  pointed  out  in  Chapter  3,  the  friction  force  may 
stop  the  relative  sliding  between  the  pair  elements  of  the 
higher  pair.  To  ensure  this  would  not  happen.  Equation  (3- 
45)  must  be  satisfied,  i.e., 

G7(X)  = u - |cosatan6|  > 0 


(5-37) 
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5. 2. 2. 3 Equality  constraint 

To  ensure  that  the  mechanism  will  not  have  an  in- 
finitely small  or  large  link  length,  it  is  necessary  to  set 
a fixed  link  length  for  either  the  input  link  or  the  output 
The  link  length  is  difficult  to  define  for  the  link 
that  attached  to  a higher  pair.  However,  for  this  par- 
ticular mechanism,  the  input  link  length  can  be  defined  as 
the  shortest  distance  form  the  input  join  axis  to  the  ini- 
tial contact  point  (or  the  center  of  the  sphere,  if 
preferred),  which  is  designated  as  P in  Figure  (5-2).  Since 
the  ground  link  length  is  already  set  to  be  unity,  we  may 
assign  P2  as  the  input  link  to  the  ground  link  length  ratio. 
Then,  an  equality  constraint  equation  can  be  formed,  as 

H8(X)  = | (u2  x P) ,x  u2 I - ?2  = 0 (5-38) 

5-2.3  Formulation  of  R-Cp-C  mechanism 

R-Cp-C  mechanism  is  shown  in  Figure  (5-4)  , which  can 
be  described  similarly  to  the  R-Sp-R  mechanism  except  that 
there  is  one  more  design  variable  to  define  the  axis  of  the 
cylinder  element  in  the  Cp  pair.  The  state  variables  are 
the  same  as  the  R-Sp-R  mechanism.  The  design  vector  is 

X = [0X,  0y,  0Z,  p,  px,  py,  pz]T  (5-39) 

where 

— — [Px'  Py/  Pz J 
a i = [ m@ ] a 
bi  = [M0]b 

£i  = a.i  x b^ 
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and 

c*j_  = cos((U)b-L  + sin(P)c^ 

is  the  unit  vector  of  the  axis  of  the  cylinder  element.  The 
angle  p is  measured  from  the  in-plane  axis  b to  the  inner 
cylinder  axis  c*.  The  rotation  matrix  [M0]  is  defined  the 
same  as  Equation  (5-28). 

5 . 2 . 3 . 1 Ob j ective  function 

There  are  four  deviation  angles  for  the  R-Cp-C 
mechanisms  as  illustrated  in  Chapter  3.  If  the  R-Cp-C 
mechanism  is  desired  to  be  a crank-rocker  type  mechanism, 
then  the  two  deviation  angles  6F23  and  6m23  will  be  a right 
angle  during  the  entire  motion  of  the  mechanism.  Thus,  it 
is  meaningless  to  maximize  their  minimum  values.  Therefore, 
the  objective  function  will  consist  of  two  other  deviation 
angles  and  the  maximum  and  minimum  accelerations,  as 

U(X)  = - k1|cos6FR|min  - k2 | cos6ft | - k3 | cosfiMR |min  + 

k4O03  + k5a03  + k6as3  + k7as3  (5-40) 

where 

CX03  = | |®3max  / ®3minl  “ ^ I > 
a03  “ I (®3max  “ ®3min) /®3max I ' 
aS3  ~ I l^3max  / ^3minl  " 1 > 

aS3  - Ms3max  " s3min^ /s3max i > 

and  03  and  S3  are  the  angular  and  linear  acceleration  of  the 
output  link,  respectively.  The  coefficients  kif  i = [1,  6], 
sre  the  weighting  factors.  The  deviation  angles 
and  6mr  are  defined  in  Chapter  3. 
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5. 2. 3. 2 Inequality  constraints 

Equations  (14)  to  (16)  are  applied  to  determine  a 
crank-rocker  type  mechanism.  The  constraint  for  continuous 
motion  is  the  same  as  the  R-Sp-R  mechanism.  Thus,  G7(X)  can 
be  similar  to  Equation  (5-37). 

5. 2. 3. 3 Equality  constraint 

Similarly  to  the  R-Sp-R  mechanism,  it  is  necessary  to 
set  a fixed  link  length  ratio  for  the  input  link  and  the 
qround  link.  The  input  link  length  is  defined  as  the  length 
of  the  common  normal  between  u3  snd  c.  Thus,  the  equality 
constraint  can  be  written  as 

H8(X)  = (P  - E>2 ) • c x u2  / |c  x u2 1 - P2  = 0 (5-41) 

where  P2  is  a given  link  length  ratio  for  the  input  link. 

5.2.4  Formulation  of  R-PT-C-R  Mechanism 

An  R-Pl-C-R  mechanism  is  shown  in  Figure  (5-6),  where 
the  input  joint  axis  U2  is  in  the  direction  of  the  x axis 
and  the  output  join  axis  u^  has  the  common  normal  (0,  1,  0)T 
with  u.2  and  has  a skew  angle  x.  A coordinate  system  a,  b, 
and  c is  attached  at  the  initial  contact  point  of  the  outer 
plane  element  of  the  PL  pair.  The  axis  a is  normal  to  the 
plane  element.  The  axis  c*  is  the  in  the  inner  plane  ele- 
ment. The  axis  of  the  floating  C joint  is  denoted  as  u3  and 
the  joint  has  a distance  vector  D3  away  form  the  input  joint 
which  is  located  at  the  origin  of  the  fixed  coordinate  sys- 


tem. 
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5. 2. 4.1  Objective  Function 


There  are  three  deviation  angles  for  force  transmis- 


sion, three  deviation  angles  for  moment  transmission,  and 


two  deviation  angles  for  the  force-produced  moment  transmis- 
sion. Thus,  the  objective  function  may  be  formed 

U(X)  = - kx | cos6f23 Iroin  ~ k2 I c°s6F34 I " k3 |cos6FR|min 
k4|cos6M23l  " k5lcos6MM34lmin  ~ k6 I cos6MRlmin 
k7  | cos6fm34  I min  “ k8  I cosSFMRlmin  + k9ct©3 
+ klO°03  + kllaS3  + k12aS3  + k13a©4  + k14a04 

(5-42) 


where 


- [0Px'  0Py'  0Pz > 0u3x'  0u3y ' 0u3z'  3'  PX'  Py/  PZ' 

*^3x'  ^3y'  ^3z^ 

= I l03max  / 03minl  “ 1| 
a03  = Me3max  “ 03min) /03max I 
aS3  ~ lls3max  / s3minl  “ 1 I 
aS3  “ Hs3max  " s3min^ /s3max  I 
a04  - I l04max  / 04minl  " 1| 
a04  = M04max  “ 04min)  /04max  I 
— i = t ^©p  ] 

b j_  = [ Mqp  ] b 

Cf  = §i  x b.j_ 

c*i  = cosfp)^  + sinfpJc-L 
^3i  = tM0u3^3 

The  coefficients  kj_,  i = [1,  14],  are  the  weighting  factors. 
[M©p]  is  the  rotation  matrix  for  the  in-plane  coordinate 
system  and  Qpj , j = x,  y,  z,  are  the  rotation  angles  for  the 
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in-plane  coordinate  system  at  the  i-th  iteration  measured 
from  the  positive  branch  of  the  fixed  coordinate  axes, 
respectively.  Similarly,  [MqU3]  is  the  rotation  matrix  for 
—3  and  ©U3  j , j = x,  y,  z,  are  the  rotation  angles  for  u.3  at 
the  i-th  iteration,  respectively.  (3  is  the  angle  measured 
from  b to  c . ©3  and  S3  are  the  relative  angular  and  linear 
accelerations  of  the  floating  C joint.  0^  is  the  angular 
acceleration  of  the  output  link.  Pj , j = x,  y,  z,  are  com- 
ponents of  the  vector  P from  the  origin  to  the  origin  of  the 
initial  moving  coordinate  system  a,  b,  c.  d3j,  j = x,  y,  z, 
are  the  components  of  the  vector  D3  from  the  origin  to  the 
floating  C joint. 

5. 2. 4. 2 Inequality  constraints 

Equations  (20)  to  (22)  are  applied  to  obtain  a crank- 
crank  type  mechanism.  Also,  Equation  (5-37)  applies  for  the 
continuous  motion. 

5. 2. 4. 3 Equality  constraints 

Since  the  R-P^-C-R  is  a four-link  mechanism,  it  is 
therefore  necessary  to  set  two  link  length  ratio  instead  of 
one  as  for  the  previous  mechanisms.  The  input  link  length 
ratio  can  be  defined  to  be  the  same  as  for  the  previous 
mechanisms.  The  second  link  length  ratio  can  be  set  between 
the  output  link  and  the  ground  link,  where  the  link  length 
is  the  length  of  the  common  normal  between  the  two  joint 
axes.  Thus,  we  have  H8(X)  to  be  defined  by  Equation  (5-37) 
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and  Hg(X)  to  be  defined,  as 

H9(X)  = (D3  - D4)  • U3  x u4  / |u3  x u4|  - P4  = 0 (5-43) 
where  P4  is  a given  link  length  ratio  for  the  output  link. 


5.2.5  Feasible  mechanism 

Feasible  mechanism"  here  means  that  the  mechanism  has 
the  ability  to  move  as  excepted,  namely,  as  a crank-crank, 
crank-rocker,  or  other  Grashof-Type  mechanisms.  The  design 
problems  formulated  are  parametric,  i.e.,  the  original  ob- 
jective function  and  the  inequality  constraint  equation  for 
the  continuous  motion  need  to  be  calculated  throughout  the 
entire  range  of  motion.  It  is  always  difficult  to  start  the 
optimization  search  from  an  infeasible  starting  point  and 
expect  that  the  constraints  will  be  satisfied  when  the 
original  objective  function  converges.  This  is  usually 
overcome  by  using  the  exterior  penalty  function.  Therefore, 
it  is  desired  to  start  the  optimization  search  from  a 
feasible  starting  point.  In  order  to  find  the  feasible 
starting  point,  an  extra  optimization  search  needs  to  be 
formulated.  This  can  be  done  by  using  the  equality  con- 
straint equation  as  the  original  objective  function.  The 
six  nonparametric  inequality  constraint  equations  for  clas- 
sifying the  type  of  mechanism  are  applied  as  the  inequality 
constraints.  Thus,  the  "pseudo-objective"  function  can  be 
written 


U 


M “ 


9 

£ Hk(X) 
k=8 


+ M (R,  Gi (X) ) 


3 = 


[1,  6] 


(5-44) 
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The  global  minimum  of  Equation  (5-44)  is  zero  when  all  the 
inequality  constraints  are  satisfied.  Then,  the  design  vec- 
tor resulting  from  this  optimization  process  will  provide  a 
feasible  mechanism,  i.e.,  a feasible  starting  point  for  the 
objective  functions  derived  in  the  previous  sections.  It 
can  be  observed  that  Equation  (5-44)  contains  non-parametric 
constraints  and  the  original  objective  function  is  much 
simpler.  Thus,  it  will  be  much  easier  to  converge  from  an 
infeasible  region  to  a feasible  region  than  by  using  the 
parametric  constraint  function  directly.  However,  the  con- 
straint for  the  continuous  motion  is  not  included  in  Equa- 
tion (5-44),  which  is  a parametric  constraint. 


CHAPTER  6 

NUMERICAL  EXAMPLES 
6 . 1 R-Sp-R  Mechanisms 

Figure  (5-1)  illustrates  the  parameters  of  the  R-Sp-R 
mechanism.  The  formulation  of  objective  function  and  the 
constraints  are  illustrated  in  Chapter  5.  The  R-Sp-R 
mechanism  is  desired  to  be  a crank-crank  type  mechanism. 

The  input  joint  is  assumed  rotating  with  unit  angular 
velocity  (©2=1  rad/sec)  and  unit  input  torque  (M^  = 1 
lbf*in).  The  static  friction  coefficient  is  set  at  y.  = 

0.27.  The  given  variables  are: 

u2  = (1.000000,  0.000000,  0.000000)T 

d2  = (0.000000,  0.000000,  0.000000)T 

u3  = (0.866025,  0.000000,  0.500000)T 

D3  = (0.000000,  1.000000,  0.000000)T 

The  above  vectors  are  illustrated  in  Section  5.2.2.  The  in- 
put link  length  is  designated  to  be  equal  to  1.5.  The  fol- 
lowing set  of  the  weighting  factors  gives  a better  over  all 
result  considering  the  displacements  and  the  deviation 
angles,  which  is 

kl  - k2  = k3  = i 

The  local  optimum  point  based  on  Equation  (5-29)  with  the 
above  weighting  factors  for  one  of  the  branches  of  the 
mechanism  is  found  at  (see  Equation  (5-27)): 
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X = [ 2.893710,  1.826581,  0.946411, 

-0.025338,  1.426154,  -0.464859]T 

The  corresponding  design  parameters  (see  Figure  (5-2))  are 
obtained  carried  out  as 

a = ( 0.247194,  -0.856662,  -0.452798)T 

b = (-0.758864,  0.119415,  -0.640207)T 

c = ( 0.602512,  0.501867,  -0.620571)T 

P = (-0.025338,  1.426154,  -0.464859)T 

The  actual  input  link  length(the  absolute  value  of  P)  is 
1.500003.  The  fixed  cosine  of  deviation  angle  (see  Equation 
(5-9))  for  the  output  force  to  produce  the  output  rotation 
is  cos6fr  = 0.958777.  The  minimum  cosine  of  the  variable 
deviation  angle  is  cos6F23  = 0*632942.  The  maximum  allow- 
able static  friction  coefficient  to  mechanism  movability  is 
0.478643.  The  minimum  plane  element  length  and  width  are 
1.511937  along  the  direction  of  the  moving  axis  b and 
1.962780  along  the  direction  of  the  moving  axis  c. 

Figure  (6-1)  shows  the  output  motions (in  radian)  of  the 
mechanism.  03 , 03 , and  03  are  the  angular  displacement,  an- 
gular velocity,  and  angular  acceleration,  respectively. 
Figure  (6-2)  shows  the  relative  displacement  of  the  Sp  pair. 
The  path  is  plotted  by  ec(the  relative  movement  along  the 
direction  of  moving  axis  c)  relative  to  e^fthe  relative 
movement  along  the  direction  of  the  moving  axis  b) . Figure 
(6  3)  illustrates  the  thansmissivity  and  the  corresponding 
deviation  angles.  The  curve  of  force-to-torque  transmis- 
sivity  the  absolute  value  of  the  scalar  product  of 
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©3=  angular  displacement  (rad) 

©3:  angular  velocity  (rad/sec) 

©3:  angular  acceleration  (rad/sec2) 


Figure  6 1 Output  motion  of  R-Sp-R  mechanism. 


ec(in) 
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Figure  6-2  Relative  path  of  the  elements  of  the 


Sp  pair. 
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Figure  6-3  Force-to-Torque  Transmissivity  and  the 
corresponding  deviation  angles. 
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c°s6F23  and  cos6fr  (see  Equation  (3-18)).  Figure  (6-4) 
shows  the  magnitudes  of  the  normal  force  and  sliding  force 
with  and  without  friction.  Figure  (6-5)  shows  the  output 
torque  with  (|Mr|)  and  without  (|Mr*|)  friction.  Figure 
(6-6)  shows  the  power  loss  which  is  the  scalar  product  of 
the  magnitude  of  sliding  relative  velocity  and  the  magnitude 
of  the  friction  force.  The  instantaneous  efficiency  (the 
ratio  between  the  output  power  and  the  input  power)  is  shown 
in  Figure  (6-7).  Since  the  input  power  is  unity,  the  in- 
stantaneous  efficiency  is  numerically  equal  to  the  output 
power . 


6 . 2 R-Cp-C  Mechanisms 

An  R-Cp-C  mechanism  is  illustrated  in  Figure  (5-4), 
which  is  designed  as  a crank-rocker  mechanism.  The  angular 
velocity  of  the  input  joint  is  unity.  The  output  joint  is 
attached  to  a translational  and  a rotational  springs.  The 
spring  indexes  for  the  translational  and  rotational  springs 
are  kt  = 0.1  lbf/in  and  kr  = 1 lbf/rad,  respectively.  The 
output  force  and  torque  are  proportional  to  the  transla- 
tional and  rotational  displacements  of  the  output  link,  as 
seen  in  Equations  (4-15)  and  (4-16).  Thus,  the  output  link 
may  contribute  either  positive  or  negative  work.  The  input 
torque  is  an  unknown  quantity,  which  may  also  be  doing  posi- 
tive or  negative  work  to  the  input  link.  The  static  fric- 
tion coefficient  is  set  at  u = 0.27.  The  given  variables 
are  designated  as  follows: 
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Fn:  normal  force  with  friction 

jlf 

Fn  • normal  force  without  friction 
fs:  sliding  force  with  friction 

Fs  : sliding  force  without  friction 


Figure  6 4 Normal  forces  and  sliding  forces  at  the  Sp  pair 
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Mr:  output  torque  with  friction 

Mr*:  output  torque  without  friciton 


Figure  6-5  Output  torques. 


129 


P.L.  : 

power  loss  (lbf* in/sec) 

Vs  = 

relative 

sliding  velocity  (in/sec 

Ff  : 

friction 

force  (lbf) 

Figure  6-6  Power  loss. 
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Figure  6-7  Instantaneous  efficiency. 
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u2  = (1.000000,  0.000000,  0.000000)T 
d2  = (0.000000,  0.000000,  0.000000)T 

u3  = (0.866025,  0.000000,  0.500000)T 
D3  = (0.000000,  1.000000,  0.000000)T 
The  input  link  length  (the  length  of  the  common  normal  be- 
tween u2  and  c)  is  set  at  0.5.  The  following  weighting  fac- 
tors set  gives  a better  overall  result  considering  the  dis- 
placements and  the  deviation  angles,  which  is 
kl  = k2  = k3  = k4  = k6  = 1,  k5  = k7  = 5 
The  local  optimum  point  (see  Equation  5-39)  based  on  Equa- 
tion (5-40)  with  the  above  weighting  factors  for  one  of  the 
branches  of  the  mechanism  is  found  at 

X = [-0. 23213805 E+00,  0 . 20730233E+01 , 0 . 32297996E+01 , 

0 . 1463 1425E+01 , 0 . 24620417E+00 , 0 . 64163337E+00 , 
0.57991137E+00]t 

The  corresponding  parameters  then  are  obtained  as 
a = (-0.036556,  -0.700632,  -0.712585)T 
b = ( 0.136410,  0.702889,  -0.698097)T 

c = ( 0.998903,  -0.046490,  -0.005534)T 
P = ( 0.246204,  0.641633,  0.579911)T 

The  input  link  link  length  is  equal  to  0.500002.  The  mini- 
mum plane  element  length  is  2.408281  along  the  direction  of 
the  moving  axis  b.  The  minimum  cylinder  element  length  is 
3.576647  along  the  direction  of  the  moving  axis  c*.  The 
fixed  cosines  of  deviation  angles  for  the  output  force  to 
produce  output  translation  and  rotation  are  cos6FT  = 

0.387951  and  cos6FR  = 


0.662558,  respectively. 
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When  the  frictional  effects  are  taken  into  account,  the 
width  of  the  plane  element  becomes  essential.  If  the  width 
of  the  plane  element  is  small,  i.e. , the  length  of  the  con- 
tact line  is  short,  the  normal  forces  produced  by  the  input 
torque  at  the  two  contact  points  at  the  edges  of  the  plane 
element  will  be  very  large  in  order  to  produce  enough  trans- 
mitting moment.  In  this  case,  the  corresponding  friction 
forces  become  too  large  such  that  the  mechanism  may  not  be 
able  to  move  at  all.  This  can  be  observed  from  the  virtual 
work  equation  such  as  Equation  (4-23).  The  power  loss  is 
proportional  to  the  absolute  magnitudes  of  the  normal  force. 
The  theoretical  power  loss  may  be  larger  than  the  input 
power  if  the  absolute  magnitudes  of  the  normal  force  are  too 
large.  Then  there  is  no  way  to  produce  positive  output 
power  in  order  to  overcome  the  output  resistant  force  and 
torque.  Also,  when  solving  the  static  equilibrium  equations 
(4-19)  and  (4-20),  there  is  no  way  to  get  positive  values  of 
the  magnitudes  of  the  reaction  force  no  matter  how  the  signs 
of  the  normal  forces  change  due  to  a moment  change.  So,  the 
width  of  the  plane  element  is  increased  until  Equations  (4- 
19)  and  (4-20)  can  yield  positive  values  for  the  magnitudes 
of  the  reaction  forces  with  proper  signs  of  the  normal  force 
vectors.  The  minimum  plane  width  is  found  at  W = 10  when 
the  diameter  of  the  cylinder  element  is  neglected.  Note 
that  if  the  diameter  of  the  cylinder  element  is  counted,  the 
width  of  the  plane  element  needs  to  be  increased.  This  is 
because  there  will  be  one  more  friction  moment  component  as 
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illustrated  in  Equation  (2-11).  Thus,  the  power  loss  will 
be  increased.  Then,  a longer  contact  line  is  needed  to 
reduce  the  absolute  magnitudes  of  the  normal  forces. 

Figure  (6-8)  and  (6-9)  illustrate  the  translational  and 
rotational  output  motions  of  the  R-Cp-C  mechanism.  03 , 03 , 
and  03  are  the  angular  displacement,  angular  velocity,  and 
angular  acceleration,  respectively.  S3,  S3,  and  S3  are  the 
linear  displacement,  linear  velocity,  linear  acceleration, 
respectively.  Figure  (6-10)  shows  the  relative  displacement 
of  the  Cp  pair.  The  path  is  plotted  with  the  quantities  ec 
and  ej-,  which  are  the  relative  movements  of  the  initial  mid- 
point of  the  contact  line.  Figure  (6-11)  illustrate  the 
f orce-to-torque  and  f orce-to-force  transmissivities  (see 
Equations  (3-18)  and  (3-21)  ) and  the  corresponding  deviation 
angles.  These  deviation  angles  are  measured  at  the  midpoint 
of  the  contact  line.  The  deviation  angles  measured  at 
thetwo  contact  points  are  illustrated  in  Figure  (6-12). 
oos6f231  and  cos6F232  denote  the  deviation  angles  for  inter- 
nal force  transmission.  cos6FF^  and  cos6FF2  denote  the 
deviation  angles  for  producing  output  rotational  motion. 
Figure  (6-13)  shows  the  deviation  angles  for  moment  trans- 
missions and  the  torque-to-torque  transmissivity ( see  Equa- 
tion (3-23)).  Figures  (6-14)  and  (6-15)  show  the  normal 
forces  with  and  without  friction,  respectively.  The  posi- 
tive and  negative  values  indicate  the  direction  of  the  nor- 
mal force.  Fnl  and  Fn2  denote  the  normal  forces  at  the  two 
contact  points  with  friction.  Fn-]_*  Fn2*  denote  the  normal 
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s3:  displacement  (in) 

i 

S3:  velocity  (in/sec) 

S3:  acceleration  (in/sec2) 


Figure  6-8  Output  translational  motion  of  R-Cp-C 
mechanism. 
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©3 : displacement  (rad) 

©3:  velocity  (rad/sec) 

• • 

©3:  acceleration  (rad/sec^) 


Figure  6-9  Output  rotational  motion  of  R 


-Cp-C  mechanism. 
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Figure 


10  Relative  displacement  of  the  initial  contact 
point ( initially  the  midpoint  of  the  contact 
line)  at  the  Cp  pair. 
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Ficfuxe;  6-11  Force-to-Force  and  Force-to-Torgue 

Transmissivities  and  the  corresponding 
deviation  angles. 
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Figure  6 12  Deviation  angles  at  the  contact  points. 
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figure  6 13  Torgue-to-torque  transmissivity  and  the 
corresponding  deviation  angles. 


140 


Figure  6-14  Normal  forces  at  the  Cp  pair  with  firctional 
effects . 
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Figure  6-15 


Normal  forces  at  the  Cp  pair  without  frictional 
effects . 
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forces  at  the  two  contact  points  without  friction.  Fn  and 
Fn  denote  the  resultant  normal  forces  with  and  without 
friction,  respectively.  In  Figure  (6-16),  |Mn|  and  |Mn*| 
are  the  magnitudes  of  the  normal  moments  at  the  midpoint  of 
the  contact  line  of  the  Cp  pair  with  and  without  friction, 
respectively.  |Mf|  is  the  magnitude  of  the  friction  moment. 

(6-17)  illustrates  the  input  power  and  power  loss. 

The  positive  and  negative  values  of  the  input  and  output 
power  represent  a positive  and  negative  work  contributed  by 
the  input  joint  or  the  output  joint,  respectively.  Figure 
(6-18)  shows  the  instantaneous  efficiency.  The  first  por- 
tion of  the  curve  is  the  ratio  of  the  output  power  at  joint 
3 to  the  input  power  at  joint  2.  The  second  portion  of  the 
curve  is  the  reverse  ratio  since  now  joint  3 becomes  the  in- 
put and  joint  2 becomes  the  output.  The  instantaneous  ef- 
ficiencies for  output  rotation  and  translation  are  not  il- 
lustrated.  The  mechanical  advantages  are  shown  in  Figure 
(6-19).  The  input  force  and  output  force  are  assumed  ap- 
plied at  unity  distance  from  joint  2 and  joint  3,  respec- 
tively. Then,  the  mechanical  advantages  are  obtained  from 
the  ratios  of  output  force  and  torque  to  input  torque, 
respectively. 


6 • 3 R-Pr.-C-R  Mechanisms 

An  r-pl-c-r  mechanism  is  shown  in  Figure  (5-6),  which 
is  designed  as  a crank-crank  type  mechanism.  The  input 
joint  is  rotating  with  unit  angular  velocity  (©2=1 
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Figure  6-16  Moments  at  the  Cp  pair. 


sec 
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Figure  6-17  Input  power(I.P.),  output  power(O.P.),  and 
power  loss ( P . L . ) . 
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Input  Angles 


Figure  6-18  Instantaneous  efficiency. 
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mechanical  advantage  for  translatioal  motion 
M.A.j-:  mechanical  advantage  for  rotational  motion 


Figure  6-19  Mechanical  advantages. 
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rad/ sec)  and  unit  input  torque  (Md  = 1 lbf*in).  The  static 
friction  coefficient  is  set  at  y.  = 0.27.  The  given  vari- 
ables are: 

u2  = (1.000000,  0.000000,  0.000000)T 

d2  = (0.000000,  0.000000,  0.000000)T 

u4  = (0.866025,  0.000000,  0.500000)T 

d4  = (0.000000,  1.000000,  0.000000)T 

The  input  and  output  link  lengths  are  set  at  0.8  and  0.4, 
respectively.  The  wighting  factors  are  chosen  as 

ki  = 1;  i=l,  2,  3,  4,  5,  6,  7,  8,  9,  11,  13 

kj  = 3;  j = 10,  12 

The  optimal  design  variables  based  on  Equation  (5-42)  with 
the  above  weighting  factors  for  one  of  the  branches  are 
found  at 

X=  [ 0 . 64428810E+00 , 0 . 22138480E+01 , 0 . 14650180E+01 

0. 28783130 E+01,  0 . 11716420E+01 , 0 . 19570740E+01 

0. 16722840 E+01,  0 . 99362500E+00 , -0 . 40688370E+00 

0 • 75471 430E+00 , -0 . 88621710E-02 , -0 . 38934150E+00 
0.88036170E+00]t. 

The  corresponding  parameters  are  obtained  as 
a = (-0.082442,  0.311164,  -0.946774)T 

b = ( 0.481539,  -0.819312,  -0.311203)T 
c = (-0.916835,  -0.396079,  -0.050339)T 
P = ( 0.993625,  -0.406884,  0.754714)T 

u3  = ( 0.046998,  0.965289,  -0.256921)T 

D3  = (-0.008862,  -0.389342,  0.880362)T 

The  input  and  output  link  lengths  are  equal  to  0.799992  and 
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0.39995.  The  minimum  inner  plane  length  along  the  direction 
of  the  moving  axis  c*  is  1.4049.  The  minimum  outer  plane 
length  along  the  direction  of  the  moving  axis  b is  3.4710. 
The  fixed  deviation  angles  6F34  and  61423  &te  found  at 
cos6p34  = 0.841836  and  cos6]423  = 0.996596.  Similar  to  the 
R-Cp-C  mechanism,  the  widths  of  the  plane  elements  affect 
the  magnitudes  of  the  normal  forces.  The  sufficient  widths 
for  both  plane  elements  are  2.4,  for  which  Equations  (4-34) 
to  (4-36)  can  yield  positive  values  for  the  reaction  forces 
with  proper  signs  of  the  normal  force  vectors. 

Figure  (6-20)  shows  the  output  motions.  The  rotational 
and  translational  motions  of  link  3 with  respect  to  link  4, 
connected  by  the  floating  C joint  are  shown  in  Figure  (6-21) 
and  (6-22),  respectively.  The  relative  displacement  of  the 
PL  pair  is  shown  in  Figure  (6-23).  The  path  is  the  relative 
movement  of  the  initial  gravity  center  of  link  2 between  the 
two  plane  elements,  refered  to  the  plane  elememt  of  link  3. 
Figure  (6-24)  shows  the  deviation  angles  for  force  transmis- 
sion at  the  three  contact  points  on  the  P^  pair  ( see  Appen- 
dix C).  Figure  (6-25)  shows  the  force-to-torque  transmis- 
sivity (see  Equation  (3-26))  which  is  the  product  of 
cos6F23'  cos6fr  and  the  invariant  cos6F34.  The  curve  of 
cos6f23  is  taken  from  the  assumption  that  the  force  is 
transmitted  at  the  gravity  center  of  the  contact  surface. 
Figure  (6-26)  shows  the  torque-to-torque  transmissivity  (see 
Equation  (3-30)  and  the  corresponding  cosines  of  deviation 
angles  6^34  and  6jy[R. 


149 


©4:  displacement  (rad) 

©4:  velocity  (rad/sec) 

©4:  acceleration  (rad/sec2) 


Figure  6-20  Output  motion  of  r-pl-C-r  mechanism. 
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©3 : displacement  (rad) 

©3:  velocity  (rad/sec) 

©3:  acceleration  (rad/sec2) 


Figure  6 21  Relative  rotational  motion  of  link  3 with 
respect  to  link  4 in  the  C joint. 
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S3:  velocity  (in/sec) 

S3:  acceleration  (in/sec2) 


Figure  6-22  Relative  translational  motion  of  link  3 
with  respect  to  link  4 in  the  C joint. 
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£c (in) 


Figure  6 23  Relative  displacement  of  elements  of  the 
PL  pair:  path  of  a point  of  link  2 with 
respect  to  link  3 . 


153 


Figure  6-24  Deviation  angles  for  internal  force 

transmission  at  the  three  contact  points. 
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Figure  6-25  Force-to-torque  transmissivity  and  the 
corresponding  deviation  angles  (cos6M/1 
0.841836). 
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Figure  6-26  Torque-to-torgue  transmissivity  and  the 
corresponding  deviation  angles  (cosSmtj 
0.996596). 
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Figure  (6  27)  shows  the  (force  torque ) ^-to-torque  transmis- 
sivity (see  Equation  (3-36)  and  the  corresponding  cosines  of 
deviation  angles  6pM34'  6fmR*  Figure  (6-28)  and  (6-29)  show 
the  transmitted  normal  forces  at  the  three  contact  points 
and  the  summation  of  these  three  normal  forces  with  and 
without  frictional  effects,  respectively. 

When  there  is  no  friction,  the  summation  of  the  three 
normal  forces  is  zero.  This  indicates  there  is  no  force 
transmitted  through  the  Pp  pair.  This  can  be  explaned  by 
the  free-body  diagram  of  the  R-PL-c-R  mechanism  (Figure  4-3) 
and  the  facts: 

1.  The  C joint  is  assumed  frictionless  so  that  there  is 
no  force  transmitted  along  the  direction  of  the 
joint  axis  (see  Equation  (4-30)). 

2.  The  normal  of  the  plane  element  can  not  be  perpen- 
dicular to  the  axis  of  the  C joint  to  prevent  the 
infinite  sliding  (see  Equation  ( F-7 ) ) . 

The  moments  taken  at  the  gravity  center  of  the  PL  pair 
are  shown  in  Figure  (6-30).  |Mn|  and  |Mn*|  denote  the  mag- 
nitudes of  the  normal  moments  with  and  without  friction. 

I ^f I denotes  the  absolute  magnitude  of  the  friction  moment. 
Figure  (6-31)  shows  the  magnitudes  of  the  output  torque  with 
and  without  frictional  effects,  the  output  power,  and  the 
power  loss.  |Mr|  denotes  the  output  torque  with  friction. 

|Mr  | denote  the  output  torque  without  friction.  Since  the 

1.  The  torque  exerted  by  the  transmitted  force  in  the  PT 
pair  induced  at  the  C joint  L 
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Figure  6-^7  Force  torgue-to-torque  transmissivity  and  the 
corresponding  deviation  angles. 
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Figure  6-28  Transmitted  normal  forces  at  the  PL  pair  with 
fricitonal  effects. 
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Figure  6-29  Transmitted  normal  forces  at  the  PT  pair 
without  f ricitonal  effects. 
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L Pair. 


Figure  6 30  Moments  at  the  gravity  center  of  the  P 
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Figure  6-31  Output  torque _ with  friction  (Mr),  output  torque 
without  friction  (Mr  ),  output  power(O.P.),  and 
power  loss ( P . L . ) . 
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input  torque  is  unity,  the  mechanical  advantage  is  numeri- 
cally equal  to  the  output  torque.  Also,  the  input  power  is 
unity,  thus,  the  instantaneous  efficiency  is  numerically 
equal  to  the  output  power. 


CHAPTER  7 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

7 . 1 Conclusions 

The  transmission  characteristics  of  spatial  mechanisms 
containing  one  higher  pair  has  been  illustrated  by  the  use 
of  deviation  angles  and  transmissivities.  Similar  defini- 
tions are  applied  to  spatial  mechanisms  containing  lower 
pairs  only.  This  study  provides  a full  understanding  of  the 
way  to  improve  transmission  quality.  The  deviation  angles 
for  internal  force  and  moment  transmissions  provide  a clear 
view  of  the  characteristics  of  each  specific  pair.  There- 
ter , the  transmissivities  link  the  corresponding  deviation 
angles  to  show  the  contributions  of  the  input  force  or 
torque  to  the  resistant  force  and/or  torque.  Thus,  the 
transmission  quality  of  mechanisms  can  be  improved  by  using 
the  criteria  developed  in  this  work  in  the  design  optimiza- 
tion. 

The  force  analysis  of  the  higher  pairs  and  some  lower 
pairs  in  Chapter  2 provides  a basic  understanding  of  these 
pairs.  The  way  how  the  forces  and  torques  are  transmitted 
is  related  to  the  geometric  constraints  of  the  pairs.  It 
has  been  shown  that  the  force-closed  higher  pairs  possess  a 
transmission  characteristic  that  needs  to  be  paid  more  at- 
tentions to  prevent  separation  of  the  pair  elements.  After 
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the  transmission  of  force  and  torque  at  each  pair  are 
analyzed,  the  deviation  angles  and  transmissivities  can  be 
defined  as  shown  in  Chapter  3.  With  all  those  criteria  and 
the  numerous  design  variables,  the  search  of  a mechanism 
with  better  transmission  quality  has  not  been  easy.  Thus,  a 
feasible  mechanism  must  be  found  before  all  the  constraints 
are  taken  into  account.  Then,  the  mechanisms  can  be  further 
optimized  according  to  the  objectives.  However,  the  ex- 
amples are  quite  general  in  their  design  purposes.  These 
examples  are  only  used  to  show  the  characteristics  of  these 
Particular  mechanisms . The  emphasis  should  be  on  the  poten- 
tial usages  of  these  design  criteria.  They  can  be  included 
as  constraints  or  part  of  the  objective  function  for  the 
design  of  other  types  of  mechanisms. 

Friction  always  plays  a significant  role  in  mechanism 
design.  When  friction  is  involved,  the  way  of  modeling  the 
reaction  forces  and  moments  becomes  important.  There  is 
more  than  one  way  to  model  the  forces  and  moments  acting  in 
specific  pairs  for  static  force  and  torque  analysis.  The 
forces  and  moments  in  this  work  have  been  modeled  in  such  a 
way  that  the  frictional  effects  can  be  obtained  directly. 
From  the  examples  in  Chapter  6 , we  observed  that  friction 
does  affect  the  dimensions  of  the  higher  pairs.  Of  course, 
friction  also  affects  the  efficiency  of  the  mechanisms.  It 
can  be  foreseen  that  with  the  consideration  of  friction  in 
each  joint,  friction  needs  much  more  careful  attentions.  If 
3-11  joints  are  considered  with  friction,  the  force  and 


165 


torques  will  have  to  be  solved  for  iteratively.  This  is 
especially  true  when  link  masses  are  included.  If  the 
masses  are  considered,  elastic  analysis  will  be  needed  for 
those  mechanisms  that  are  statically  indeterminate. 

The  main  objective  of  this  work  is  to  look  into  the 
transmission  characteristics  and  it  has  been  done  by  the  use 
of  deviation  angles  and  transmissivities.  Including  the 
static  force  and  torque  analysis  exposes  the  significant 
frictional  effects.  It  is  hoped  that  this  work  fulfilled 
the  absence  of  the  definitions  of  transmission  criteria  in 
spatial  mechanisms  containing  higher  pairs.  It  is  also 
hoped  that  this  work  brought  more  attentions  to  the  area  of 
frictional  effects  in  mechanism  design. 

7 . 2 Recommendations 

With  the  design  criteria  of  transmission  quality  having 
been  defined,  including  these  criteria  in  the  synthesis  pro- 
cedures becomes  possible.  These  synthesis  procedures  can 
include  motion  generation  or  path  generation. 

The  examples  in  this  work  are  limited  to  mechanisms 
whose  configuration  has  two  branches  only.  Otherwise,  it 
would  be  difficult  to  optimize  the  mechanisms  to  satisfy  the 
constraints  during  the  entire  range  of  motion.  For  those 
mechanisms  having  more  than  two  branches,  the  avoidance  of 
branching  and  sequence  defects  becomes  essential. 

As  mentioned  in  the  previous  section,  dynamic  analysis 
including  the  frictional  effects  is  a worthwhile  objective 
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for  further  research.  However,  it  can  be  expected  that  a 
certain  degree  of  difficulties  will  be  encountered. 


APPENDIX  A 

RELATIVE  DISPLACEMENTS  OF  AN  Sp  PAIR 

The  contour  of  the  plane  in  the  Sp  pair  has  been 
derived  by  Hernandez  (1983),  which  also  includes  the  design 
of  the  inner  sphere  and  its'  inserted  rod.  In  this  work, 
the  detailed  design  of  the  pair  elements  is  not  considered 
except  for  the  physical  size  of  the  plane  element.  It  is 
necessary  to  limit  the  size  of  the  plane  element  in  the  op- 
timization  process  since  the  contour  of  the  plane  may  be 
very  large  to  satisfy  some  of  the  constraints.  Therefore,  a 
procedure  similar  to  that  illustrated  by  Hernandez  to  deter- 
the  contour  of  the  plane  element  is  adopted  here.  A 
coordinate  system  consisting  of  a,  b,  and  c axes  is  attached 
to  the  plane  where  the  origin  is  located  at  the  initial  con- 
tact point  P,  as  shown  in  Figure  (A-l).  Let  a be  the  normal 
of  the  plane  element;  b and  c^  are  on  the  contacting  plane. 
The  j-th  position  vectors  of  the  initial  contact  point 
derived  by  starting  out  from  the  input  and  output  joint  axes 
are  denoted  as  Pj  and  Pj ■ , respectively.  Then  we  may  write 
the  following  equation 

£j  = £j'  + ebjbj-  + ecj£j'  (A-l) 

where  e^j  and  ecj  denote  the  relative  displacements  between 
the  contact  point  along  b j i and  c j ■ axes,  respectively;  here 
— j ' and  £j ' are  obtained  by  starting  out  from  the  output 
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A 

b 

I 


Figure  A-l  Relative  displacements  of  the  Sp  pair 
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joint  axis,  since  the  plane  is  assumed  to  be  connected  to 
the  output  joint.  The  scalars  e^j  and  ecj  may  be  obtained 
by  taking  the  dot  product  of  Equation  (A-l)  with  b j . and 
Cj . , respectively,  which  gives 

ebj  = (£j  “ Pj  ' ) * bj  . (A-2) 

and 

ecj  = (Pj  - Pj  « ) • Cj  . ( A-3  ) 

The  difference  of  the  maximum  positive  and  negative  values 
of  ebj  gives  the  extent  along  axis  b of  the  plane  element 
needed.  Similarly,  the  extent  along  axis  c is  determined  by 
the  maximum  positive  and  negative  values  of  ec j . The  plot 
of  the  two  scalars  e^j  and  ecj  gives  the  contour  of  the 
plane  element. 


APPENDIX  B 

RELATIVE  DISPLACEMENTS  OF  A Cp  PAIR 

The  line  of  action  of  the  force  in  the  Cp  pair  can  be 
considered  to  pass  through  the  axis  of  the  cylinder.  Thus, 
the  contact  line  can  be  assumed  to  be  parallel  to  the  axis 
of  the  cylinder.  Figure  (B-l)  shows  the  initial  and  j-th 
position  of  the  mid-point  P of  contact  line  c*.  A mooving 
coordinate  system  consists  of  the  normal  of  the  plane,  a, 
and  the  vectors  b and  c.  The  vector  c*  is  the  axis  of  the 
cylinder  element  and  forms  an  angle  3 with  b.  Thus,  we  may 
consider  c to  be  rotated  about  a by  the  relative  rotation 
angle  3 formed  with  b.  Then  we  may  write 

c*  = cos^b  + sinf3(a  x b)  (B-l) 

The  two  scalars  ecj  and  e^j  define  the  required  extent 
of  the  plane  and  length  of  the  cylinder  to  prevent  separa- 
tion at  this  j-th  position.  Then,  we  may  have  the  relation- 
ship of  the  position  vectors  of  point  P expressed  as 

£j'  + ebjkj'  = P j + e*cj£*j  ( B-2 ) 

where  for  the  case  of  P2-CP-P2  three-link  mechanism  with  the 
plane  connected  to  the  input  link  and  the  cylinder  connected 
to  the  output  link,  and: 

£*j  = 2*  + (cos02  - l)U2c*  + sin©2 (u2  x £* ) (B-3) 

bj.  = b + (cos03  - l)u3b  + sin03(u3  x b)  (B-4) 

Pj  = P + (cos02  - 1)U222  + sin02(u.2  x 0.2 ) (B-5) 
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Figure  B-l  Relative  displacements  of  the  Cp  pair. 
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P j i = P + (cos©3  - DU323  + sin03(u,3  x £3 ) (B-6) 

where 

]±2c*  = (H2  x £*)  x u2; 

U3b  = ( H3  x b)  x u3; 

Uigi  = (Hi  x 2i)  x Hi; 

2i  — ?.i  - 2i  f i = 2 , 3 

Pj  and  c*j  are  the  j-th  position  vector  of  the  contact  point 
and  the  unit  vector  of  the  direction  of  the  cylinder  element 
axis,  respectively,  derived  by  starting  out  from  the  input 
axis  (i.e.,  a point  fixed  to  the  plane  that  is  connected  to 
the  input  link).  P j 1 and  b j 1 are  the  j-th  position  vectors 
of  the  same  contact  point  and  the  unit  vector  parallel  to 
the  U-shape  element,  respectively,  derived  from  the  output 
axis  (i.e.,  a point  fixed  to  the  inner  plane  which  is  con- 
nected to  the  output  link) . P is  the  initial  position  vec- 
tor of  the  contact  point;  is  the  vector  from  the  origin 
to  the  axis  u^;  9j_  is  the  relative  rotation  of  the  i-th 
joint . 

By  taking  the  scalar  product  of  Equation  (B-2)  with 
(£  j x P j ) and  (bj 1 x P j 1 ) respectively,  gives  the  two 
scalars  as 


eb  j 

and 


x Pj 


x P- 


( B-7 ) 


e 


* 

cj 


x £ j ' 
x P j ' 


( B-8 ) 
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The  contour  of  the  plane  element  can  be  obtained  by- 

plotting  the  scalar  e^j  along  the  b axis  and  sin(3j£*cj  along 

c axis.  The  largest  values  of  e^j  and  e*cj  determine  the 
extent  of  the  plane  and  the  minimum  length  of  the  cylinder, 
respectively.  If  the  width  of  the  plane  is  denoted  as  W, 
the  end  points  of  the  contact  line  can  be  expressed  as 

£lj  = POj  + (Wj/2)c*j  (B-9) 

£2 j = ^Oj  ~ (Wj/2)c*j  ( B— 1 0 ) 

where 

P0j  = P j + e*cj£*j  = P j ' + ebjMj ' and 

Wj  = W/sin|3 j . 


APPENDIX  C 

RELATIVE  DISPLACEMENTS  OF  A PL  PAIR 

A procedure  similar  to  that  in  Appendix  B can  be  used 
to  derive  the  position  of  the  contact  points  of  the  PL  pair 
(see  Figure  C-l).  The  same  coordinate  system  is  used  to 
derive  the  two  scalars  ebj  and  e*cj  as  in  Equations  (B-3, 
B-4).  The  lengths  of  the  inner  plane  and  the  U-Shaped  outer 
plane  element  are  determined  by  the  positive  maximum  values 
of  e*cj  and  ebj , respectively.  The  position  of  the  four 
possible  contact  points,  i.e,  the  four  corners  of  the  con- 


tact surface,  can  be  expressed  as 

£lj  = P0j  + (Hj/2)bj.  + (Wj/2)c*j  (C-l) 
P.2  j = loj  + (Hj/2  )bj  i - (Wj/2)c*j  ( C-2 ) 
P3j  = P0j  * (Hj/2 )bj . - (Wj/2)c*j  ( C-3 ) 
P4  j = P0j  " ( Hj /2 )bj 1 + (Wj/2)c*j  ( C-4 ) 


where  Pg j , Wj  ,and  c*j  are  defined  the  as  in  Appendix  B;  and 
— j ' = — j3  + ( cos©4  - l)U4bj3  + sin04(u.4  x bj3) 
kj3  = k + (cos©3  - l)U3b  + sin03(u3  x b) 

Hj  = H/sin^j 

The  quantity  H is  the  width  of  the  inner  plane  element,  W is 
the  width  of  the  U-Shaped  outer  plane  element.  Then,  the 
three  contact  points  as  used  in  Chapter  2 for  modeling  the 
reaction  force  and  moment  can  be  arbitrarily  chosen  from  the 
four  corner  points.  Similar  to  the  Cp  pair,  the  contour  of 
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L Pair. 


Figure  C-l  Relative  displacements  of  the  P 


176 


the  U-shaped  outer  plane  element  can  be  plotted  by  the  two 
scalars  e^j  and  sin£je*cj  which  are  the  relative  displace- 
ments along  the  in-plane  coordinates  b and  c,  respectively. 


APPENDIX  D 

DISPLACEMENT  EQUATIONS  OF  R-Sp-R  MECHANISM 

Figure  (D-l)  shows  the  j-th  position  of  an  R-Sp-R 
mechanism.  The  constraint  equation  is 

(Pj ' " ) * a j . =0  (D-l) 

where 

P j = £ + (cos©2  - 1)U2q2  + sin02(u.2  x £2) 

P j ' = P + (cos©3  - DU323  + sin©3(u.3  x ) 

a j 1 = a + (cos©3  - l)U3a  + sin©3(u.3  x a) 

and 

^iQ3  = ^i  x 2i)  x ui;  2i  = £ " Dj.  i = 2,  3 

^3a  = (ii3  x x u3 . 

By  rewriting 

P j ' = 23  + ( cos©3  - l)U3g3  + sin©3 ( u3  x Q3 ) + D3  (D-2) 
which  can  be  denoted  as  23 j 1 , the  constraint  equation  can  be 
simplified  to  be 

23  * ^ + (P.3  " Pj  ) * aj  . = 0 (D-3) 

Then,  the  above  equation  can  be  rearranged  into  a quadratic 
equation  in  terms  of  X3 , where  X3  = tan-1(©3/2),  cos©3  = (1 
- X32)/(l  + X32),  and  sin©3  = 2X3/ (1  + X32),  as 

aX32  + bX3  + c = 0 (D-4) 

where 

a = (P  - Pj)  • a - 2(D3  - Pj)  • u3a 
b = 2 ( D3  - P j ) • u.3  x a 
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Figure  D-l  The  R-Sp-R  mechanism  in  the  j-th  position. 
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c = ( P - P j ) • a 

The  angular  velocity  ©3  is  obtained  by  taking  time 
derivative  of  Equation  (D-l)  with  ©3  constant,  which  gives 


0 


3 


where 


( D-5 ) 


P'j  = - sin02U222  + cos©2(u2  x g2) 
a'ji  = - sin©3U2a  + cos©3(u.3  x a) 

P'jt  = - sin©3U323  + cos©3(u.3  x 23 ) 

By  taking  time  derivatives  twice  of  Equation  (D-2) 
gives  the  angular  acceleration  of  the  output  link,  as 

..  - ( d©3^  + e©2©3  + f©2^) 

03  = ( D-6  ) 

g 

where 

d = P"j.  • a-j.  + 2 (P ' j » • a'j.)  + (Pj,  - P j ) • a"  j . 
e = - 2P ' j -a'j. 

f = - r j • 4j  • 

9 = £'  j • • ?_j  • + (Pj  • - Pj  ) • a'  j . 

and 

P"ji  = - cos©3U323  - sin©3 ( u3  x 23) 
a"j>  = - cos03U3a  - sin03(u3  x a) 


APPENDIX  E 

DISPLACEMENT  EQUATION  OF  R-Cp-C  MECHANISM 

Figure  ( E— 1 ) shows  the  j-th  position  of  the  R-Cp-C 
mechanism.  The  constraint  equations  are  written 

(Pj ' - Pj ) * aj .=  0 ( E-l ) 

aj . • c*j  = 0 ( E-2 ) 

where 

P j ' = P + (cos03  - DU.323  + sin03(u.3  x 23)  + S3U3 
P j = P + (cos02  - 1)U2£2  + sin02(u.2  x £2 ) 
a j 1 = a + (cos©3  - l)U3a  + sin©3(u.3  x a) 

£*j  = 2*  + (cos02  - l)U2c*  + sin02(u2  x c * ) 

and 

U3Q3  = (M3  x 23)  x M3 ; H2Q2  = (M2  x £2)  x M2? 

M3 a = (M3  x M)  x M3?  M2c*  = (M2  x £* ) x uj; 

£i=P"Mi  i = 2,  3 

By  the  use  of  the  tangent  of  the  half  angle  of  03 , 
Equation  (E-2)  can  be  written  into  a quadratic  equation. 

aX32  + bX3  + c = 0 (E-3) 

where 

a = (a  - 2U3a)  • c* j 

b = 2U3  x a • c* j 
* 

c = a • c j 

The  coefficients  a,  b,  and  c are  functions  of  the  input 
angle  ©2. 
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Figure  E-l  The  R-Cp-C  mechanism  in  the  j-th  position. 
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The  sliding  displacement  S3  can  now  be  found  from  Equa- 
tion ( E-l ) . This  gives 

(Pj  - P-5* ) • a-;  i 

S3  = — ( E-4 ) 

U3  * a j . 

where 

Pj*  = P + (cos©3  - DU323  + sin©3(u.3  x £3) 

From  the  denominator  of  Equation  (E-4),  we  can  observe  that 
U3  can  not  be  perpendicular  to  a j . . Otherwise,  the  sliding 
of  the  Cp  pair  may  go  to  infinity. 

The  angular  velocity  ©3  is  obtained  by  taking  time 
derivative  of  Equation  (E-2)  with  @2  constant,  which  gives 


where 

c*  ' j = - sin©2U2C*  + cos©2(u.2  x c * ) 
a'ji  = - sin©3U3a  + cos©3(u.3  x a) 

The  translational  velocity  S3  is  then  solved  by  taking 
time  derivative  of  Equation  (E-4)  with  ©2  constant.  This 
leads  to 


©2 ( P ' j * a j , ) + @3 [ ( P j - Pj.)  • a ' j . - P'j.* 
S3  = 

^3  * a j ' 

where 

P'j  = - sin©2.U22i2  + cos©2(u.2  x 22^ 

P'j  1*  = - sin©3U323  + 00503(113  x 23) 

The  angular  acceleration  is  given  as 


] 


( E-6 ) 
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- [032 ( a" j i • c*j)  + 2©203 ( a ' j » • c*'j)  + 022(aj.  • c*Mj) 


where 

a"ji  = - cos03U3a  - sin03(u3  x a) 
c*"j  = - cos02U.2c*  ~~  sin02(u,2  x 2*) 

The  linear  acceleration  of  the  output  link  is 

..  - [d03  + e03^  + ©3  ( f ©2  + gS3  ) + g022] 

s3  = 

—3  * a j . 

where 

d = P'j«*  • a j . + (Pj . - P j ) • a ' j i 
e = P"j,*  • a j . + 2 ( P ' j . * • a'ji)  + (Pj . - Pj 
f = - 2P'j  . a' j . 
g = u3  • a' j . 
h = - P"j  • a j . 

and 

P"ji*  = - cos03U323  - sin03(u3  x 23) 

P"j  = - cos0]_U121  “ sin0]_(u.3_  x Qi ) 

P'ji*  = - sin03U323  + cos03(u3  x 23) 

P'j  = - sin02U.222  + cos02(u.2  x 0.2 ) 


( E— 7 ) 


( E-8 ) 


APPENDIX  F 

DISPLACEMENT  EQUATIONS  OF  R-PL-C-R  MECHANISM 

The  finite  screw  displacements  for  the  three,  four,  and 
five-link  mechanisms  containing  Sp,  Sg,  Sc,  and  Cp  pairs 
have  been  derived  by  Ghosal  (1982),  which  does  not  include 
the  P^  pair.  However,  the  displacement  equations  for  the 
mechanisms  containing  a P^  pair  can  be  derived  in  a similar 
manner  as  in  Ghosal 's  work.  Figure  (F-l)  shows  an  R-PL-C-R 
mechanism  in  the  j-th  position. 

The  PL  pair  has  one  rotation  about  the  normal  sectors 
a.j  and  a.j  ' between  the  two  plane  elements  where  the  two  nor- 
mal sectors  must  be  the  same.  The  other  relative  motion  be- 
tween the  pair  elements  of  the  P^  pair  is  the  translation  on 
the  contact  plane,  which  indicates  that  the  direction  sector 
of  the  relative  motion  must  be  perpendicular  to  the  normal 
of  the  plane.  Thus,  the  constraint  equations  for  the  planar 
pair  are  expressed  as  follows: 

' (F-l) 

(Pj ' “ £j ) * a j = 0 ( F-2 ) 

where 

a j = a + (cos02  - l)U2a  + sin02(u.2  x a) 

-j ' = §.j3  + (cos04  - l)U4aj3  + sin04(u4  x a-j3) 
aj3  = a + (cos03  - l)U3a  + sin©3(u3  x a) 

P j = P + (cos©2  - DU2Q2  + sin02(u.2  x 2.2 ) 
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Figure  F-l  The  R-PL-C-R  mechanism  in  the  j-th  position. 
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-j ' “ -j3  + ( cos©4  - 1 ) — 4Q j 4 + sin©4 (U4  x Qj4) 

+ s 3~3  j 4 

Pj3  = P + (cos03  - DU323  + sin03(u,3  x 23) 

2j4  = £j3  “ 

By  taking  scalar  product  with  u4  of  Equation  (F-l), 
gives  one  equation  correlating  02  and  03  as 

[a  + (cos03  - l)U3a  + sin03 ( u3  x a)  - a j ] • u4  = 0 

( F— 3 ) 

With  the  use  of  tangent  of  the  half  angle  formulations  using 
tan (03/2)  Equation  (F-3)  can  be  rearranged  into  a quadratic 
form  as 

aX32  + bX3  + c = 0 (F-4) 

where 

a = ( §:  ~ fij  - 2U3a ) • U4 
b = 2(u3  x a)  • u4 
c = ( a - a j ) • U4 

Then,  the  output  angular  displacement  04  can  be  ob- 
tained by  taking  the  scalar  product  of  Equation  (F-l)  with 
LUaj  3 and  u4  x aj3  to  eliminate  sin03  and  (cos©4  - 1), 
respectively,  as 

cos04  = d/e  ( f— 5 ) 

sin04  = f/e  (F-6) 

where 

d = a j • ^4aj3 
e = 1 - ( a j 3 • U4 ) 2 

f = a j • u.4  x a j 3 

Then,  04  = tan  -1-  ( sin04/cos04 ) . Note  that  in  programming. 
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the  function  AT AN 2 should  be  used  in  stead  of  ATAN  to  avoid 
getting  an  improper  sign  of  ©4. 

The  linear  displacement  of  joint  3,  S3,  now  can  be  ob- 
tained directly  from  the  other  constraint-equation  (F-2),  as 
(Pj  - Pjt*)  • aj 

S3  = ( F-7  ) 

—3 j 4 * aj 

where 

£j'*  = lj3  + (cos04  - l)U42j4  + sin04(u4  x £j4) 
u.3  j 4 = M.3  + ( cos04  - l)U4u3  + sin04(u4  x u3 ) 

Similar  to  the  R-Cp-C  mechanism,  the  denominator  of  the 
above  equation  tells  that  U3j4  can  not  be  perpendicular  to 
a j i . Otherwise,  an  infinite  sliding  will  occur. 

To  derive  the  angular  velocities,  we  can  take  the  time 
derivative  of  Equation  (F-l)  with  03  constant,  which  gives 
the  following  equation 

03a3ji  + 04a4j . = 02a' j (F-8) 

where 


-3j'  ~ a ' j3  + (cos04  - l)U4a. j3  + sin04(u4  x a'j3) 
a4j  . = - sin04U4^j3  + cos04(u4  x a.  j 3 ) 


a'j3  - - sin03U3a  + cos03(u3  x a) 

« • 

Then,  03  and  04  are  given  by  taking  dot  products  of  Equation 
(F-8)  with  u3  x a4j  i and  u.3  x a3 j 1 , respectively,  as 

U3  x a4 j . )©2 


®3  = 


(a'j 


a3 j 1 * u.3  x a4 j 1 


( F-9 ) 


and 


* (a'j  • u.3  x a3  j i ) 02 

04  - 

& j 1 • U3  x j i 


( F-10 ) 
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where 

a'j  = - sin©2U.2a  + cos©2(u.2  x a) 

The  sliding  velocity  of  joint  3,  S3,  can  be  obtained  by 

taking  time  derivatives  of  Equation  (F-2),  which  gives 

02g  - [03P3  j ' + ®4  ( P4  j ' + S3u ' 3 j 4 ) ] • a-j 

S3  = — — ( F— 1 1 ) 

^3j4  * *j 

where 

g = P'  j * aj  + (Pj  - Pj . ) - a'j 
P'j  = - sin©2U222  + cos©2(u.2  x 2.2 ) 

P4j  i = - sin©4U42j4  + cos©4(u.4  x Q. j 4 ) 

P3  j ' = P'j3  + (cos©4  - 1 ) U4P ' j 3 + sin©4(u.4  x P'j3) 

P'j 3 = - sin©3U323  + cos©3(u.3  x 23) 
u'3j4  = - sin©4U4U3  + cos©(u_4  x u3 ) 

The  angular  accelerations  ©3  and  ©4  can  be  obtained  by  first 
taking  the  time  derivative  again  of  Equation  (F-8),  which 
gives 

©3a3 j ' + ©32a33ji  + 2©3©4a34j . + 04a4 j . + ©42a44j . 

= ®22^"j  (F-12) 

where 

a"j  = - cos©2U2a  - sin©2(u.2  x §) 
a"j3  = - cos©3U3a  - sin©3(u.3  x a) 

a33 j 1 = a"j3  + (cos©4  - l)U4aiij3  + sin©4(u.4  x a"j3) 

_a34j  1 = - sin©4U4aij3  + cos©4(u.4  x a. ' j 3 ) 
a44j  i = - cos©4U4aj3  - sin©4(u.4  x a j 3 ) 

Then,  by  taking  dot  products  with  u3  x a3  j ■ and  u.3  x a4j  1 , 
respectively,  gives 
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and 


03  = 


u3  x a4 j i 
• u3  x a4j i 


( F-13 ) 


04  = 


where 


u3  x a3  j i 

* U3  x a3j. 


( F-14 ) 


h = 022a"j  - 032a33jt  - 203©4a34jt  - 042a44 j . 

The  relative  linear  acceleration  in  joint  3 is  obtained  by 
taking  the  time  derivative  again  of  Equation  (F-ll),  which 
gives 

S3  = {[022P"j  - ©32P33j<  - ®42(P44j<+  S3u"3j4)  - 

"®3— 3 j ' - ©4  ( P4  j 1 + S3u'3j4)  - 2©304P34j,  - 
20403H' 3 j 4 ^ * aj+  [2©2(02P'j  " ®3l3 j ' - ©4(P4j<  + 
S3M' 3 j 4 ) " S3u3j4]  • a ' j + ©2(Pj  - Pj.*  - 
^ 3—3  j 4 ) * a"j>  / (y.3  j 4 * aj) 

where 


P"j  - - cos02U222  + sin02(u2  x 22) 

P33j,  = P"j3  + (cos04  - l)U4p"j3  + sin04(u4  x P"j3) 
P34j ' = ~ sin04U4pij3  + cos04(u4  x P ' j 3 ) 

P44j.  = - cos04U4g|j  4 - sin04(u4  x Qj4) 

L*"  3 j 4 = “ cos04U4u3  - sin©4  (u4  x u3 ) 

— " j 3 = " COS03U323  " sin©3(u3  x £3 ) 


APPENDIX  G 

ROT AT AB I L I TY  OF  R-Sp-R  MECHANISM 

The  constraint  equation  for  the  R-Sp-R  mechanism  can  be 
written,  as 

' " £ j ) • a j i =0  (G-l) 

where  P j , P-j  ■ , and  a j ■ are  denoted  in  Appendix  D.  By  the 
use  of  tangent  of  the  half  angle  of  the  output  joint  rota- 
tion angle,  ©3,  the  above  equation  can  be  rearranged  into  a 
quadratic  equation  of  X3 , where  X3  = tan(©3/2): 

aX32  + bX3  + c = 0 (G-2) 

The  coefficients  of  Equation  (G-2)  are  listed  in  Appendix  D. 
The  two  roots  of  Equation  (G-2)  correspond  to  the  two  pos- 
sible values  of  the  output  angle  03 . When  the  input  link 
reaches  a limit  position,  the  two  output  angles  become  the 
same,  i.e.,  the  determinant  of  Equation  (G-2)  must  be  equal 
to  zero.  The  discriminant 

D = b2  - 4ac  (G-3) 

can  be  rearranged  into  a quartic  equation  in  terms  of  the 
output  angle  by  the  use  of  tangent  of  the  half  angle  for- 
mulations for  the  input  rotation  angle,  which  gives 

d4x24  + d3x23  + d2x22  + dlx2  + d0  = 0 ( G-4 ) 

where 

d4  = e0  ~ 2e^_  + 4e4 
d3  = 2e2  “ 4e3 
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d2  = 2eg  - 2e]_  + 4e5 

dl  = 2e2 
d0  = e0 

and 

e0  = (23  ' x a)2 

el  = 2(23  *113  x a)(U2g2  ’ 1±3  x a)  + 

2(23  * U3a)(U222  * §) 
e2  = 2(23  *1*3  x a)(u2  x 22  ' H3  x a)  + 

2(23  ‘ H3a ) (li2  x 22  ‘ a) 

e3  = 2(U222  * H3  x a) (u2  x 22  * U3  x a)  - 2(U2£2  * 

a)(u2  x 22  ' a)  + 2(U222  * a)(u2  x 22  * U3a)  + 

(J±2  x 22  * a)(U222  • U3a) 

e4  = (^222  * H3  x a) 2 - (U2S2  * a) 2 + 

2 ( —222  * a)(U2g2  • U3a) 
e5  = (u2  x 22  * H3  x a) 2 - (u2  x 22  * a)2  + 

2(u2  x 22  * a) (u2  x 22  * U3a> 

Solving  Equation  (G-4)  will  yield  four  values  of  the 
input  angle  corresponding  to  the  four  possible  limit  posi- 
tions. However,  if  an  input  crank  link  is  desired,  no  limit 
positions  should  exist,  therefore.  Equation  (G-4)  should 
have  no  real  roots. 

Similarly,  to  ensure  a fully  rotatable  output  link,  we 
may  write  the  out-input  equation  by  reversing  the  roles  of 
the  input  and  output  rotation  angles  in  Equation  (G-2), 
which  gives 

fX22  + gx2  + h = 0 (G-5) 


where 


192 


f — §3  * (°lj  ' ~ a)  - 2U2Q2  * aj  ' 

9 = u2  x Q2  • a j . 
h = 23  ‘ (aj i - a) 

The  discriminant  of  Equation  (G-5)  can  be  rearranged  into  a 
quartic  equation  in  terms  of  X3 , as 

i4X34  + i3X33  + i2X32  + i]_X3  + i0  = 0 (G-6) 

where 

14  = 30  " 2 j 1 + 4 j 4 
i3  = 2 j 2 - 4 j 3 

j-2  = 2 j 0 ~ 2 j 1 + 4 j 5 
Li  = 2 i 2 ; io  = jo 

and 

jo  = (y.2  x 0.2  * a)2 

jl  = 2 f (y.2  x 22  * a)(u2  x 22  ‘ U3a)  + 

^—222  * a)  (23  * U3a.)  ] 
j 2 = 2 [ ( u2  x 22  * a ) ( u2  x 22  * u3  x a ) + 

^—222  * a) (23  * U3  x a) ] 
j 3 = 2t(u2  x 22  * H3a)(u2  x 22  * 113  x a)  - 
(23  ‘ L[3aM23  * U3  x a)  + (23  * U3^) 

(—222  * )±3  x a)  + (23  * U3  x a) (U2£2  * U3a) 
j4  = (u2  x 22  * U3a)2  " (23  * U3a)2  + 

2(23  • M3a)(y222  * L!.3a) 
j5  = (Hl2  x 22  * H3  x a)2  - (23  • U3  x a)2  + 

2(23  ‘ U3  x a)(U222  * H3  x a) 

If  a fully  rotatable  output  link  is  desired.  Equation  (G-6) 
must  have  no  real  roots  so  that  there  are  no  limit  positions 
of  the  output  rotation  angle. 
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